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We study the low temperature phase of the 3D Coulomb glass within a mean field approach which 
reduces the full problem to an effective single site model with a non-trivial replica structure. We 
predict a finite glass transition temperature Tc, and a glassy low temperature phase characterized by 
permanent criticality. The latter is shown to assure the saturation of the Efros-Shklovskii Coulomb 
gap in the density of states. We find this pseudogap to be universal due to a fixed point in Parisi's 
flow equations. The latter is given a physical interpretation in terms of a dynamical self-similarity of 
the system in the long time limit, shedding new light on the concept of effective temperature. From 
the low temperature solution we infer properties of the hierarchical energy landscape, which we use 
to make predictions about the master function governing the aging in relaxation experiments. 
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I. INTRODUCTION 

In Anderson insulators the Coulomb interactions be- 
tween localized electrons remain essentially unscreened 
and give rise to a strongly correlated low temperature 
state, the so-called electron or Coulomb glass. This state 
of matter is expected to occur in strongly doped, but 
insulating semiconductors, in granular metals, as well 
as in dirty thin metal films. Such systems were long 
ago predicted to exhibit glassy propertiesii^ii^ii^ due to 
their inability to reach the ground state on experimen- 
tal timescales. The latter leads to experimentally ob- 
servable out-of-equilibrium phenomena such as the slow 
relaxation^ii of conductivity and compressibilit y^'^'^° , 
agingiiii^ii^ as well as memory effects^'i^'i^. The change 
of current noise characteristics across the metal insulator 
transition has also been ascribed to glassy behavio r^^i^^ . 
Further, a series of experiments in doped samples close 
to the metal insulator transition have shown a variety of 
glassy featuresi^i^i^'^. 

In recent years there has been growing experimental 
evidence^ that the glassy behavior observed in some of 
the above disordered electronic systems is an intrinsic 
property of the interacting electrons themselves, suggest- 
ing that those systems are indeed realizations of Coulomb 
glasses. In addition there is ample numerical evidence for 
glassy behavior in such system a^'^i^^'^^i^^'^ ''i^^ . 

From a different point of view, insulators with strong 
Coulomb interactions have long been known to exhibit 
a prominent suppression in the density of states around 
the chemical potential, as first argued by PoUak^S and 
Srinivasan22,. Later, Efros and Shklovskii"^^ have shown 
that a pseudogap p(e) < Ce^^^ in the local density 
of states is a necessary prerequisite for any configu- 
ration of a D-dimensional Coulomb glass to be sta- 
ble with respect to single particle hops. The pres- 
ence of this pseudogap leads to a substantial increase 



of the resistivity from Mott's variable range hopping 
R{T) « i?oexp[(rM/T)i/(^+i)] (assuming a constant 
density of states) to Efros- Shklovskii's law, R{T) « 
i?o exp[(TEs/T)^/^]. Furthermore, it significantly sup- 
presses tunneling at low voltages, as was demonstrated 
experimentally'^^ . 

For a long time the relation between these two aspects 
of strongly interacting Anderson insulators has remained 
unclear. However, a close analogy with spin glasses sug- 
gested that there may be a deeper connection between 
theroli23i3i. Indeed, Efros' and Shklovskii's (ES) stabil- 
ity argument for the Coulomb gap can also be applied to 
long range spin glasses, in particular to the Sherrington- 
Kirkpatrick (SK) model. In that case the argument leads 
to the conclusion that " the density of states" , or more ap- 
propriately, the distribution of local fields acting on the 
spins, must at least have a pseudogap p(e) < |e| to en- 
sure stability with respect to the simultaneous fiip of two 
spins. On the other hand, it is known that for the SK 
model the presence of a linear pseudogap at low temper- 
ature is related to the occurrence of a finite temperature 
glass transition and the ensuing ergodicity breakin g^^i^^ . 
These observations hold true also for a fully connected 
electronic model (a fermionic SK model with small quan- 
tum fluctuations)^ as pointed out by Dobrosavljevic and 
Pastor. These authors further suggested that a similar 
relationship might hold between Coulomb gap and the 
"Coulomb glass" phase in systems with l/r interactions. 
Miiller and loffe"^- have recently shown that this is indeed 
true, as will be analyzed more thoroughly in the present 
paper. 

The ES stability argument with respect to single par- 
ticle hops or spin flips only provides an upper bound 
for the density of states or the distribution of local fields, 
while it is difficult to argue rigorously for the actual satu- 
ration of this minimally required power-law suppression. 
The same is a priori true for long ranged spin glasses. 
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However, in the case of the SK model, one can actu- 
ally prove the saturation of the linear bound. Bray and 
Moore found that there are massless modes in the exci- 
tation spectrum around typical minima in phase spacers, 
and that this marginal stability almost implies the pres- 
ence of a linear pseudogap in p(e). The latter was finally 
directly shown to be a property of the exact low temper- 
ature solution of the SK model^i^. 

For the case of electron glasses, a recently introduced 
replica mean field approach^^^S?, can be used to describe 
Coulomb correlations in a self-consistent local approxi- 
mation. This approach not only predicts a finite tem- 
perature glass transition, but similarly as for the SK 
model, an analysis of the low temperature solution shows 
that the marginality of the glass phase guarantees the 
presence of a saturated Efros-Shklovskii Coulomb gap 
p{E) (X E^^^ at low energies^. 

In this paper, we present a detailed study of this mean 
field theory for 3D electron glasses and discuss its pre- 
dictions for experiments and numerical studies. We will 
frequently compare the formalism for the electron glass 
to analogous concepts established in spin glasses, and in 
particular for the SK model. After the introduction of 
the model and a review of the mean field approxima- 
tion in Sec. |TT1 we discuss the high temperature phase in 
Sec, mil and introduce the notion of thermodynamic and 
instantaneous local densities of states which will be at 
the focus of the low temperature study. 

In Sec. IIVI we show that a glass transition occurs due 
to strong fiuctuations in the Thomas-Fermi screening 
which destroy the exponential screening present in the 
high temperature phase. This effect leads to a critical 
(marginally stable) and correlated glassy state through- 
out the whole low temperature phase. In this phase, the 
replica symmetry is spontaneously broken, signaling the 
presence of many metastable states. The near degener- 
acy of the latter is at the origin of the marginality whose 
physical implications will be discussed. 

In Sec. |V] we study the glass phase in detail. We give a 
physical interpretation of the formalism required to solve 
Parisi's replica symmetry breaking scheme in the low 
temperature phase, and prove the permanent marginal- 
ity of the glass. In Sec. IVII we obtain the local density 
of states as a function of temperature and disorder from 
the self-consistent mean field solution. In particular we 
will see how the phase transition is related to the opening 
of the Coulomb gap below Tc- Sec. I VIII focuses on the 
asymptotic low temperature regime and establishes that 
the local density of states exhibits the Efros-Shklovskii 
pseudogap with a saturated exponent D — 1. Further, 
we discuss how lattice effects modify this asymptotic low 
energy result at intermediate energies to make the gap 
exponent appear larger (« 2.34 instead of 2 in 3D). 

In Sec. IVIIII we review the generalization of the 
fluctuation-dissipation theorem and the occurrence of 
an effective temperature, and show how it is encoded 
in Parisi's ultrametric Ansatz and the replica formal- 
ism. From the presence of an asymptotic fixed point 



in the Parisi's renormalization group-like fiow equations 
we obtain a new local interpretation of the effective tem- 
perature, and infer a self-similar structure of the long 
time dynamics. The latter is found to be accompanied 
by a gradual decrease of the average screening length, 
rsc{t) - e^/KT,tf{t). 

As we show in Sec. IIXI the structure of the low tem- 
perature solution yields non-trivial information about 
the structure of phase space in the vicinity of a given 
metastable state. Combining it with a trap model for 
dynamics we discuss the possible relevance to aging ex- 
periments. Further possible manifestations of the glassy 
phase in experiments are addressed in the Discussion 
(Sec|X|. We discuss limitations of mean field theory and 
compare the 3D results with numerical approaches, sug- 
gesting various future studies to test our predictions. The 
main results are summarized in the Conclusion fSec. lXl| l. 



II. SINGLE SITE MODEL FOR THE COULOMB 
GLASS 

A. Disorder average of the lattice Coulomb glass 

We consider the classical lattice model for Coulomb 
glasses, 

where rii G {0, 1} describes the occupation of site i, K is 
the average site-occupation, ei is a random on-site poten- 
tial and K is the host dielectric constant. We restrict our 
analysis to the case of a half-filled cubic lattice, K = 1/2, 
and a Gaussian distribution of the disorder potential. 
Both assumptions are not crucial to our results, but sim- 
plify the further analysis. As long as disorder is much 
stronger than the interaction between nearest neighbors, 
the restriction to a cubic lattice should also not be impor- 
tant since the effective low-energy sites will be randomly 
placed. However, as we will see in Sec. I VIII lattice ef- 
fects persist down to relatively low energies, and thus 
very large disorder is necessary for the glass transition to 
become insensitive to the underlying lattice. 

For the following, it will be convenient to use the Ising 
spin notation Si — 2ni — 1 e {±1}- We consider a cubic 
lattice with spacing £ = 1 which fixes the unit of length. 
Replicating the system n times, and averaging over the 
disorder, we obtain the replica Hamiltonian 

-. n n 

a— 1 i^j a, 6—1 i 

with 

J^, = ^^ W=W^"\ (3) 
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where the overbar indicates the disorder average. Here 
we have chosen e'^/4:K£ = 1 as the unit of energy. T 
denotes a n x n matrix with all entries equal to 1. In 
order to describe a quenched disorder average, the num- 
ber of replicas has to be sent to n ^ eventually. The 
replicated spin problem ([2|) is amenable to standard high 
temperature expansions in the Coulomb interactions, as 
outlined in App. [X] 



B. Mapping to a single site problem 

In the limit of strong disorder, one finds that the ir- 
reducible vertex insertions of the high temperature ex- 
pansion of App. [A] are dominated by site-diagonal terms 
while off-diagonal contributions are suppressed by higher 
powers oil/W. This suggests to look for a self-consistent 
mapping to an effective single site model which allows us 
to resum the family of leading diagrams with local vertex 
insertions i^^i^ 

From a physical point of view large disorder suppresses 
the Thomas Fermi screening by quenching the electrons. 
This preserves the long range of the Coulomb interac- 
tions and hence a large number of effective neighbors 
for a given site. This situation is favorable for a cav- 
ity or mean-field approach where the behavior of a single 
site and its environment are described in a self-consistent 
manner. 

On a technical level, a similar approximation to the 
above "locator" approximation is made in dynamical 
mean field theories^, which is usually justified by in- 
voking the limit of large dimensionality. In the present 
problem, however, we can use the long range nature of 
the essentially unscreened Coulomb interaction to jus- 
tify the mean field approximation. Since we are treating 
a classical, thermodynamic rather than a quantum dy- 
namic problem, we do not have a frequency dependence 
of the replica diagonal part of the self-energy. Instead, 
we obtain a non-trivial structure in the replica sector. It 
is well-known that the latter encodes information about 
dynamic behavior in the asymptotic long time limit, as- 
suming a generic Langevin dynamics of the spins^. 



C. Self-consistent description 

We seek to map the full lattice problem onto a self- 
consistent single site model with the effective Hamilto- 
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nians^ 



/3i^0({Sa}) ^-^Y. ^-(^"fc + W^^ab)sb, (4) 



whose exact solution is supposed to resum all local 
(cactus-like) diagrams of the original lattice problem. In 
order to make the mapping self-consistent we require that 
both the local irreducible polarizability 11^6 (i.e., the ir- 
reducible two-leg vertex of the diagrammatic expansion. 



cf. Eq. (jA12p ). and the local two-point functions (si^aSi^b) 
(i.e., the overlap Qab in spin terminology) be the same 
in the two models. As derived in detail in App. [A] this 
leads to the self-consistency conditions 



Qab — {SaSb)Ha 
N ^ 



(5) 



i,ai>i,b/H„ — 



(27r)3 



where brackets denote a thermal average and J'k is the 
Fourier transform of the Coulomb interaction ([3]), J7fc = 
47r/fc2 for k ^ 0. Eq. (O holds for aU (a, 6), with the 
obvious constraint for the diagonal elements Qaa = Q ^ 
1. 

Using the first line in ([5]) to substitute (3/Tl by Q^^ + 
/3^A, we obtain the self-consistency equation relating 
overlap Q and effective coupling A, 



Qab 



{SaSbjHo 



d-^k 

(2^ 



We rewrite Eq. ^ in more compact form, 
/?Q = Xi[(/3Q)"'+/3A], 
introducing the response function 

d^k 1 



Xi(p^) 



(2^)3 Jfc+p2- 



(6) 



(7) 



(8) 



Notice that the argument takes the place of an in- 
verse polarizability operator. In the context of Coulomb 
interactions it has the interpretation of the square of an 
average Thomas- Fermi screening length, whence the sug- 
gestive notation p. 

For a strongly suppressed polarizability, p ^ 1, e.g., 
due to large disorder W or at low temperature, one ob- 
tains the asymptotic behavior 



Xi(p') 



1 



p^ 



Sk 1 Ji 



(27r)3 p4 + p2 

2-^7? Clatt , 
p3 p4 



(9) 



We have used the fact that the absence of self-interactions 
ensures J{r = 0) = J^, Ju = 0. On a lattice the integrals 
over k are cut off on short scales fc ^ in such a way 
as to preserve the normalization J^, = J dP'k/{2nY — 
1. Note that the second term in Q is dominated by 
the long range properties at small k of the interactions, 
Jk 'iir/k'^ and is thus universal, i.e., independent of the 
considered lattice. The higher order terms are sensitive 
to short length scales, the leading corrections being given 
by 



latt 



fceBZ 



(fik 
(2^ 



(2^ 



= lim 



fceBZ, fc>fco 



(fk 
(2^ 



(10) 
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where BZ indicates the first Brillouin zone of the consid- 
ered lattice. Below we will consider the case of a cubic 
lattice for which one finds Ciatt = —8.9555 using the lat- 
tice Fourier transform of App. |B] 



III. HIGH TEMPERATURE PHASE 

At high temperatures, the system is in a unique ther- 
modynamic state (a single ergodic compartment), as re- 
flected by the unbroken replica symmetry of the overlap 
and coupling matrices in the self-consistent local approx- 
imation {Qa^b = Qrs, Ao^^f, = Ars). 

The physical content of this high temperature solution 
is rather easy to understand: Due to the random onsite 
fields, Ei, each site carries an average charge or "magneti- 
zation" in spin language, qi = (si)/2 = mi/2. Therefore, 
even at high temperature, two equilibrated copies a and 
f3 of the same disordered sample have a non-vanishing 
average similarity, i.e., a positive overlap Q°^^ of their 
site occupation pattern. 



\ " a 



(11) 



In the single site approximation this translates into a 
finite replica off-diagonal overlap Q"^^ — Qa^b- 

The Coulomb repulsion from the average charges qi 
augments the original Gaussian disorder on other sites 
and makes the distribution of local fields wider. The 
self-consistent approximation treats this extra disorder 
as a Gaussian distribution with width Ars • This physics 
is contained in Eqs. (j6l7p which can be rewritten for the 
replica symmetric case in the form 

Qrs = (saSb)Ho = y d2/-PRs(2/; ARs)tanh^(/3?/),(12) 
Ars = V [{rTf] Ors, (13) 

where we have introduced the distribution of effective 
local fields, Prs(2/), and the screening radius, r^^, 

p I . ^ exphj;V2(^2 + ARs)] 

^'Rs(y;ARs = ======= — , (14 

V27r(M^2 + Ars) 

{^^f = xrM/3(l~QRs)]. (15) 
We have also defined the function 

1 1 



where 



X2(p^)^-X'i(p2) 



(2^)3 (Jfc+p2)2- 



(16) 



(17) 



■0 contains the crucial information about the Coulomb 
interactions on the original lattice. A physical interpre- 
tation of this function will be given further below. 



The asymptotics of ipip^) for large screening radius p 
is obtained as 



(V) 



1/2 



2Ci 



latt 



o(l). 



(18) 



The exponent fi — 1/2 of the leading term reflects the 
universal long range tail of the Coulomb interactions in 
3D. For general long range repulsive interactions J^{r) oc 
in D dimensions it generalizes to fi — l — v/{D — v). 
Note that the high temperature solution of Eqs. (|12ll3p 
determines self-consistently the Gaussian extra disorder 
Ars. More insight into the meaning of (fT5|) can be 
gained from the analogous equation for the high tem- 
perature phase of the Sherrington-Kirkpatrick spin glass 
in a random field which we outline in App. [C] For that 
model one finds the exact identity '0(p^) = !> and conse- 
quently the overlap and effective coupling matrices coin- 
cide, Ars — Qrs. This is physically transparent consid- 
ering the explicit expression of the variance of the addi- 
tional random fields created by the site magnetizations 
Too, 



A 



RS 



TO,- 



^RS- 



(19) 



The case of the electron glass is more involved be- 
cause the equivalent expression with bare interactions 
Jij = l/j^jj diverges. Instead, one has to take into ac- 
count that interactions are screened by other charges, 
however, without double counting direct interactions and 
screening contributions. The locator approximation pro- 
vides a compact way to solve this problem in a self- 
consistent manner. Indeed, a formal expansion of the 
self-consistent Ars in powers of J shows that the low- 
est order term A^g — f^{JijmjJikm}~) is modified by 
screening corrections 



(^Ars = ^{Jijmj)Jikl3{l - rnl)Jki 



(20) 



which have to be resummed to yield a finite result in the 
self-consistent local approximation. 



A. The density of states 

The above high temperature formalism yields a Gaus- 
sian distribution Prs(j/) of effective local fields y. We re- 
fer to them as "thermodynamic" fields in the sense that 
the average occupation of a site is given by (si) = rrii = 
tanh(/3?/i). This field distribution controls the compress- 
ibility (charge susceptibility) 



Kc = 13(1 - Qbs) ^13 dy 



Pnsjy) 
' cosh^ (Py) 



(21) 



The field yi describes the energy to fiip the occupation 
on the site i, including local rearrangements of neigh- 
boring particles which tend to lower the energy cost as 
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compared to an isolated flip. Note that in a glassy state 
such local relaxations differ from a global relaxation for 
which large portions of the system need to rearrange 
to best accommodate the flipped site. More precisely, 
one should think of yi as taking into account relaxations 
within the phase space valley corresponding to the cur- 
rent local metastable state. 



B. Onsager term 

For certain experiments in strongly insulating mate- 
rials, such as photoemission or tunneling from a broad 
junction, what matters is rather the distribution of 
"instantaneous" local fields hi = dH/drii = ei + 
e'^ /KTijUj (the cost of isolated flips). While multi- 
particle processes may in principle occur, the cost in 
quantum action associated with the hopping of electrons 
other than the one which is emitted or tunnels largely 
exceeds the corresponding gain in density of states. The 
contribution of multi-particle processes to such observ- 
ables is thus greatly suppressed in strongly insulating 
classical Coulomb glasses. 

With some additional reasoning, we can obtain the dis- 
tribution of "instantaneous" local fields from the above 
formalism, too. In spin language, the time-averaged lo- 
cal field (hi) — — J2j=£i 'JijiTT-j: and the "thermodynamic" 
field (iyi = tanh~^(TOj) differ in general, because the lat- 
ter includes local relaxation processes in the environment, 
while the first describes solely the average energy cost to 
flip the occupation of the single site i. Consequently, if 
we hold the spin so fixed at site 0, the average local field 
it sees is larger than the thermodynamic field yo by the 
polarization response of the environment. 



(ho) 



Vo + soho- 



(22) 



The term ho is known as Onsager's back reaction, 
familiar from similar equations obtained by Thouless- 
Anderson-Palmcr (TAP) for the SK-model'^^. In that 
model, ho is given by the sum over all local cactus 
diagrams starting with two legs at site 0. Those are 
simply resummed by hp = J2i^o '^oi^i^ with the irre- 



ducible polarizability 11^ w drrii/dhi = j3{\ 



(up 



to terms of order 1/A^). Using J^^ = 1/N and the 
Edwards Anderson self-overlap Qea = ^^^X^i'^f; one 
finds ho = /?(! - Qea)- 

In the case of Coulomb glasses, the expression for ho 
has to be generalized since the equivalent of the above 
expression diverges for J'ij — 1/rij. Indeed, we have to 
sum over all higher order ring diagrams-*^, 

ho = J^iji^jiJ^jii — '^ijl^jl'^jlj2^j2'^j2i + ■ ■ ■ 

(23) 

where the IIj's are local (site-dependent) irreducible po- 
larizabilities in the presence of the fixed spin at site i. 

In the spirit of the present mean field approach, we 
approximate IIj by their site average which is justified by 



the large number of sites contributing to the reaction. To 
leading order in large disorder and at low temperatures 
we have Tlj ~ l/''sc- We may then perform the sum over 
ring diagrams and obtain 



Tr- 



20F 



J 



(24) 



with the Thomas-Fermi screening radius rgc, as intro- 
duced in Eq. (|15p for the high temperature phase. It 
will be generalized later to the low temperature regime 
(cf., Eq. (|60p ). The approximation is valid for a strongly 
suppressed susceptibility, where Tsc is large. In App.lPjwe 
derive the more precise expression ho = /3{Aaa—^b^a) — 
Xi^i'^c) — 1/kC) taking properly into account the exclu- 
sion of site i from the ring diagrams. Note that for the SK 
model where A = Q, this form reduces to the expression 
given above. 

In the high temperature, replica symmetric phase of a 
strongly disordered sample the irreducible polarizability 
is approximately equal to the bare local susceptibility 
n w xioc ~ and thus 



R.S 



2(27r) 



1/4 



(25) 



When the Onsager back reaction becomes comparable to 
the temperature, the spins develop strong correlations 
with their environment. We will see in the next section 
that these correlations actually induce a glass transition 
at Tf^ h^^ . 



C. Distribution of instantaneous fields 

With Eq. we obtained the distribution of thermo- 
dynamic fields, P{y), from the self-consistent solution of 
the mean field approach. In this subsection, we show how 
to obtain the distribution of local instantaneous fields hi, 
using the insight from the previous subsection. For the 
joint distribution of hi and the orientation of the spin s^, 
we have by definition 



1 ^ f dX 



iXh — iX J^ijS 

3 

iXh - (^j5j)s,= 



X exp 



'^V i^j^k) cg.^s '^f'i 



3,k 



(26) 



We have only retained the first two cumulants, to be 
consistent within the locator approximation. From the 
TAP-like equations we may identify the first cumu- 
lant as 



{J^ijSj) 



V^3/si=s 



(^i>s.=s = y» + she 



(27) 



The second cumulant is almost insensitive to the value of 
the spin at site i, and evaluates to 



iTij i^j^k) c g.^g Jki 

j-k 



J- 



J 

hp 

T 



1 



J 



(28) 



(29) 



in the locator approximation. A closer analysis of the 
Onsager term (cf. App. [D|, shows that this relation is 
essentially exact. 

The righthand side of Eq. can now be expressed in 
terms of averages over the joint probability distribution 
for {y,, s,), given by P(y) exp(s^?;)/2 cosh(/3?/), 



P{h,s) 



dyP{y) 



e^p{sf3y) 
2cosh(/3y/2) 



(30) 



dX 
2^' 



^i\{h-y-sho)^~\^ho/2l3 



dyP{y) 



exp[ 



p{h-y-shoY 
2ho 



2 cosh(/32/) ^2TTho/l3 



Summing over s = ±1 we finally obtain the desired dis- 
tribution of instantaneous fields, 



P{h) = / dyP{y) 



cosh(/3ft,) exp[ 



I3{h-y) 
2ho 



Phg_] 
2 i 



cosh(/3y) ^2TTholP 



(31) 



This relation generalizes an expression derived for the 
SK model in Ref. H. Notice that for fields h,y:^ ho,T 
the two distributions essentially differ by a shift, P{h) ^ 
P{y^h-ho). 



D. Correlations in the instantaneous field 
distribution P{h) 

In the high temperature phase the distribution of ther- 
modynamic fields P{y) is a featureless Gaussian, cf., 
Eq. (fT4|l . Yet, due to the Onsager back reaction, a cor- 
relation hole starts developing in the distribution of in- 
stantaneous fields P{h), well before the transition to a 
strongly correlated glassy state occurs at Tc- The high 
temperature correlation hole in P{h) is thus not directly 
related to the low temperature glass phase, but simply 
reflects particle-hole correlations in the liquid Coulomb 
"plasma". This was discussed in detail in Ref. 39, which 
we follow here in referring to this correlation hole as a 
"plasma dip". Note that while at the transition noth- 
ing particular happens to P{h), Tc marks the opening 
of a pseudogap in P{y) which is related to the breaking 
of ergodicity, as we will see in more detail in the next 
section. 

The high temperature distribution P{h) peaks at finite 
fields of the order of ho (see Fig. [T|). The peak at neg- 
ative fields is due to sites i occupied by electrons which 




FIG. 1: Evolution of the correlation hole (plasma dip) in the 
distribution of the instantaneous fields P{h) for moderate dis- 
order W — 2. The curves are plotted for the inverse tempera- 
tures /3 = 0.25,0.5,1,2,4,7 (from top to bottom at ft = 0) 
above the glass transition {/3c — 7.149). The plasma dip 
starts opening as a perturbation of order 0{Tc/T) at high 
temperature and develops into a significant depletion already 
for T > Tc 



attract holes to neighboring sites, inducing in turn a pos- 
itive potential on the site i. This is again a manifestation 
of Onsager's back reaction. The same reasoning applies 
to the peak at positive fields with the role of holes and 
electrons interchanged. These correlations are present 
at all temperatures, but they are smeared by thermal 
fluctuations when T ^ . Note that the apparent 
parabolic shape of the plasma dip above the glass tran- 
sition has little to do with the universal Efros-Shklovskii 
gap p{E) ~ which we will see emerging in the low 
temperature glassy phase, cf. Fig. [H 

The plasma dip persists also in strong disorder. In fact, 
the replica symmetric mean field theory predicts a uni- 
versal shape for P{h) in the limit of very strong disorder 
and temperatures of the order of T > Tc- This prediction 
relies on the suppression of screening in the presence of 
large disorder, and applies to instantaneous fields in the 
regime h ^ Tc, much below the scale of nearest neighbor 
interactions. Using Eq. ([3T|) and anticipating = 3Tc 
(cf. Eq. (|57)) ). vahd for large disorder, one finds 



(32) 



277^/3 



cosh(7 h) 
cosh(7 y) 



2 6 ^ 



The universal shape of P{h) is illustrated in Fig. ^ 
where we plot the scaled function P{h) = WP{h = 
h/VW), evaluated at the glass transition 1/Tc = Pc = 
^VW (with 7 given in Eq. ([55)1 below). 
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hW 



FIG. 2: Universal scaling of the plasma dip in the distribution 
of instantaneous fields in large disorder at T = Tc. 



IV. GLASS TRANSITION 

The high temperature phase of the electron liquid ex- 
hibits a glass instability when the temperature becomes 
comparable to the Onsager term ho- At large disorder, 
one finds the scaling Tc ^ ho 1/VW. This is identi- 
cal to the estimatc'^^ of the width of the Efros-Shklovskii 
Coulomb gap, which follows from equating the suppos- 
edly saturated pseudogap p{e) ^ to the bare density 
of states^ po ~ 1/W. This illustrates the fact that both 
phenomena are based on the strong electron-electron cor- 
relations building up below the scale Tc. 

An analogous estimate of the critical temperature for 
the SK-model in strong random fields predicts a transi- 
tion at Tc ~ ho ~ l/W^j which is indeed confirmed by 
the exact expression for the Almeida-Thouless instability 
(see App. [C]). 



A. The glass instability 

The instability at Tc can be understood in various 
ways. Technically, the simplest way to find the insta- 
bility is to study the free energy of the effective single 
site model. 



./3F(A) = i^[ln(-/32A-/3/n)]^^ 

a 

-^E[/ ln(/3Jfc-/?/n) 



(33) 



-In 



where 11 = 11 (A) is determined from the self-consistency 
Eq. jU, equivalent to dF/dll = 0. This free energy ex- 
pression can be obtained either by " integrating" the sad- 
dle point equations (O, or from a Baym-Kadanoff func- 



tional, restricted to diagrams with a purely local self- 
energy. 

One finds that the RS solution becomes unstable to 
replicon fluctuations around the replica symmetric high 
temperature solution (Aab — A^^+SAab with J2b ^^ab — 
and 6Aaa = 0) when 



P(y;ARS) 
dy — -r— = 1. 



' cosh** (Py) 
In the large disorder limit one obtains 

Pc = jVw, with 7 ^ 3/2(27r)i/'' 
Upon restoring lattice units this reads 



Tc - 



V2 " 6(2/7r)i/4 ' 



(34) 



(35) 



(36) 



where vq — (214^^'^)^^ is the "bare" density of electronic 
states Ei. Using Eqs. (|6l7p it is straightforward to estab- 
lish that in the limit of large disorder 



Tc = /igVS for W ^Qo, 



(37) 



conflrming the previous assertion that the glass transition 
temperature scales like the high temperature Onsager re- 
action. 



B. Critical fluctuations and breakdown of 
homogeneous screening 

A more physical understanding of the instability can 
be obtained from an analysis of the connected spin-spin 
correlation function {siSj)c of the original lattice model 
before disorder averaging. According to the above ar- 
guments, its high temperature series is dominated by 
chain-like screening diagrams with local polarizability in- 
sertions Hi, 



(38) 



^•^ij + E Jil^lJlj — 



n. 



In large disorder, we have approximately 11^ « [3{\ — 
{si)^) « P/ cosh^(/3ei), where q is the local disorder of 
site i, renormalized by the extra disorder due to Hartree 
interactions with other sites (described by tree-like dress- 
ings of the bare vertices) . To leading order in large disor- 
der, the renormalized disorder is essentially the same as 
the thermodynamic field yi, with distribution Prs (|14p . 

Calculating the disorder average of the geometric se- 
ries (j38l) yields an exponential decrease of the average 
correlations, with distance, 

— — - — -— — „ ,—2 f expfifcr 
(s(O)s(r)), = p-^n I - 



p-^n 



■2 e 



-r/B.o 



n 



(39) 



8 



where 11 « 2/-\/27rW^ denotes the polarizabiHty aver- 
aged over the (renormaUzed) local disorder e. The av- 
erage screening radius is given by Rq = I/AttH ~ W 
which grows with disorder. The self-consistent approxi- 
mation in the high temperature phase resums diagrams 
renormalizing the local polarizabilities and predicts Rq = 
(r^^)^/47r, which improves on the subleading corrections 
in 1/W. 

The above result may be misleading, however, since the 
correlation functions can fluctuate strongly and might 
not even have a definite sign. In order to capture such 
effects one should compute the variance of the correlation 
function. This requires the summation of the geometric 
series of double-lined chain diagrams depicted in Fig. O 
where the lines stand for disorder averaged single propa- 
gators 



n 



and the local insertions are given by g = II^ — 11 

1 



(40) 



(41) 



where G2{k) — J (Pr exp[ikr]Gi{r). Unlike the sum in 
Eq. ([38]) . this series becomes critical at a finite tempera- 
ture Tc where 



gG2(fc = 0) = l. 



(42) 



We can check that in the large disorder limit this coin- 
cides with the critical temperature determined from the 
instability of the single site model. The variance of the 
polarizability is dominated by 



Pnsiy) 
cosh'*(/3y) 



(43) 



Further, we have 
G2(fc = 0) 



q2 



k (1 + JkW 



= ^(1 + 0(75)), 
vn 



(44) 



in an expansion for small 11 . To this order (which probes 
the universal properties of the long range Coulomb in- 
teractions) this expansion coincides with the expansion 



for i/)[(; 



'0[l/n] (cf., Eq. [T5)l. while a more careful 



analysis shows that the two expressions actually agree to 
higher order than suggested by the leading corrections in 
(|18l44p . This shows the physical equivalence of the insta- 
bility of the glass transition ([M]) and the onset of critical 
fluctuations in the screening, Eq. (|42)) . As a byproduct, 
we obtain an interpretation of the function as the long- 
distance limit of the ladder of doubled propagator lines, 
which describes the propagation of spin fluctuations. 



< SjSj > 




+ 



FIG. 3: (Color online) Typical ladder diagrams contributing 
to the variance of the spin-spin correlation function. The sin- 
gle lines represent disorder averaged propagators, the vertices 
correspond to the variance of the irreducible polarizability. 



The criticality of the geometric series for (s(O)s(r))^ 



(s(O)s(r)) translates into a power law decay of the 
square of the correlations with distance at the critical 
point. This critical behavior is the direct analog of the 
divergence of the non-linear susceptibility in long range 
spin glasses. In the SK-model with random fields, the 
dominant part of the non-linear susceptibility is given by 
the sum ^"^X^ij {^i^j)'c- K diverges as the Almeida- 
Thouless instability is approached due to the same kind 
of criticality as discussed above. 



C. Nature of the glass transition 

The Landau expansion for the single site model in 
terms of fluctuations 5K has the same structure as that 
of the random held SK-model. For both models the 
glass transition is continuous in the sense that the effec- 
tive coupling matrix A, as well as physical properties like 
the density of states, evolve continuously across the tran- 
sition. Furthermore, the breaking of ergodicity is weak 
around Tj., the conflgurational entropy (the logarithm of 
the number of metastable states) increasing continuously 
from zero at Tc- This is in contrast to structural glasses 
which usually undergo a discontinuous jamming transi- 
tion. 



V. THE GLASSY PHASE 
A. Parisi's Ansatz 

Below Tc the replica symmetry of the single site model 
is broken spontaneously which reflects the occurrence of 
many metastable states and the associated dynamical 
slowing down. 

In order to describe the low temperature phase, we take 
Parisi's ultrametric Ansatz for the pattern of the replica 
symmetry breaking in the glassy phase^i^. Technically, 
this amounts to grouping the n ^ replicas into sets 
of < xi < 1, and X2/X1 {x2 < xi) such sets into a 
bigger clusters of X2 replicas etc., building an ultrametric 
tree structure among all replicas. The values assigned 
to the coupling matrix A^b are taken to be a decreasing 
function of the distance on this tree. The coupling matrix 
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may thus be conveniently parametrized as a piecewise 
constant and increasing function A(a::) with < a; < 1. 
In the present case, the optimal self-consistent solution 
will turn out to tend to an infinite number of hierarchical 
subdivisions (continuous replica symmetry breaking) as 
described by a continuous coupling function A{x), as in 
the SK model. Accordingly the overlap between replicas 
translates into a monotonously increasing function Q{x). 



Parisi's overlap function (cf., Ref. '53' for a recent discus- 
sion of this point). 

The dynamics described by Sompolinsky's formalism 
can be understood as a slow relaxation which starts from 
an out-of-equilibrium configuration which is in general 
arbitrarily far away from the ground state. This in- 
terpretation is consistent with the picture arising in the 
mean field aging dynamics as obtained by Cugliandolo 
and KurchaUi ^^i^^ 



B. Static and dynamic interpretation 

Parisi gave an equilibrium interpretation of this 
Ansatz^, showing that P{q) = dx / dQ\Q(^x)=q describes 
the probability of two configurations of the system, sam- 
pled with the Boltzmann weight, to have mutual overlap 
q. This interpretation has recently been generalized be- 
yond the mean field case^. 

An alternative interpretation of Parisi's Ansatz which 
will prove useful for the following, was given by Sompolin- 
sky^ in terms of out-of-equilibrium Langevin dynamics. 
He supposed the existence a hierarchical set of time scales 
^1 <C ^2 *C ... ^ tfc, all of which are much larger than the 
microscopic time scale to. On the level of mean field 
equations for the dynamics, he assumed an anomaly in 
the response at each time scale, modifying the standard 
response predicted by the fluctuation-response theorem 
by a multiplicative factor x(ti) < 1, 



T 



(45) 



where R{t) is the integrated response function and Q{t) 
is the local spin-spin correlator averaged over time t and 
over all sites. The expression TdR{t)/dx{t) is sometimes 
referred to as the response anomaly A{t). Later work^ 
has clarified the meaning of this relationship as arising 
from sampling different metastable states in response to 
an external field, which is sensitive to the configurational 
entropy (number of states), rather than the standard en- 
tropy, entailing the factor x. We will come back to this 
point in Sees. IVIII Bl and IIXI where we will discuss the 
violations of the fluctuation dissipation theorem and ag- 
ing. 

However, Sompolinsky's theory was criticized because 
it invoked the divergence of all time scales and their ra- 
tios ti+i/ti with the system size N, which is inconsis- 
tent^. The theory makes sense, however, if we think 
of the time scales U as diverging, e.g., with the wait- 
ing time tuj after a sudden quench, rather than with N. 
In the limit where a strong hierarchy of times exists, one 
can show that the only time dependence of observables is 
captured by the value of the correlations, or equivalently, 
by the decreasing function x{t), which thus can serve as 
a parametrization of time. Remarkably, in the limit of an 
infinite number of time scales fc — > oo the average spin- 
spin correlation function, (si(t)sj(O)) = Q{x(t)) = Q{x), 
of the dynamical formalism turns out to be identical to 



C. Solving the self-consistent equations below 

In this subsection we show how the self-consistency 
equations ([S][7]) can be solved under the assumption that 
the replica matrices Aa^, Qah assume a hierarchical Parisi 
form. 

In order to compute {saSbjHo in ® we follow the 
scheme introduced by Sommers and Dupont"^^. We in- 
troduce a free energy per replica 4>{x, y) corresponding 
to the Sompolinsky time scale x (implicitly defined as 
the time scale over which the system explores metastable 
states up to an overlap distance Q{x) from the initial 
state), and in the presence of an external field y, 

exp[a;/30(a;, y)] = ^ exp [^(2:, y, {sq})] , (46) 

n{x,y,{Sa}) = Y X! ^''•^^ab^b + P^ySa- 
a, 5—1 a— 1 

Here PS^jj denotes the x x a;-matrix Kab — A(a;) with 
a,6G {l,...,a;}). The representation ([46]) allows for the 
derivation of a recursion equation for (j)(x—dx, y) in terms 
of 4>{x,y) (see, e.g., Ref. ISTi ). In the limit of continuous 
A(a;) and dx 0, they reduce to Parisi's differential 
equation 



dx 



IdK 
2~dx 



dy 



2 



dy 



(47) 



The boundary condition at x = 1 follows from (j46p as 

</)(x = 1, y) = A - A(l) + - log [2 cosh(/3y)] . (48) 

P 

It is often more convenient to work with the derivative 
m{x,y) — d(j3{x,y)/dy, which satisfies the flow equation 



dm IdK f d'^m ^ ^ d[inn? 



-(3x- 



(49) 



dx 2 dx \ dy'^ dy ) ' 

and initial conditions on shortest time scales 

•m{x — l,y) = tanh(/3y). (50) 

The above is a priori a formal algebraic trick 
to compute the free energy density (f)(0, 0) = 
lim„^o exp[/377o] — l]/n, from which observables can 
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be obtained by derivatives. However, it can be given a 
deeper physical meaning in the spirit of Sompohnsky's 
dynamical interpretation. As shown by Sommers and 
Dupont"^^, m{x,y) can be interpreted as the mean occu- 
pation (magnetization in spin language) of a site, aver- 
aged over the "time scale" x and in the presence of a 
local field y which is frozen over that time scale. 



D. Distribution of local fields 

We further introduce the distribution of these time- 
averaged local fields, produced by interactions with 
other parts of the system and remaining frozen on the 
"timescale" x due to the glassiness. This distribution 
P{x,y) is implicitly defined by the requirement that 



for all ai, . . . , flr G {1, . . . , x}. In the continuous limit, 
the ensuing recursion relations relating x and x + dx lead 
to the flow equations 



dP 

dx 



IdA 
'2~dx 



dy"^ 



2(3x 



d[P m] 
dy 



(51) 



The boundary condition for the largest time scales {x 
0) is 



P(x = 0,y) 



exp{-j/V2[PF^-}-A(0)]} 
^27r[T^2 + A(0)] ' 



(52) 



which describes the self-consistent Gaussian disorder pro- 
duced by the bare disorder as well as by Hartree inter- 
actions with the frozen part of charges that persists to 
the longest time scales. As in the replica symmetric case 
A(0) does not vanish because of the random field disor- 
der. The same holds true in the random field SK model. 

The object of ultimate interest is the distribution of 
frozen local fields P{x = l,y), the low temperature ana- 
log of the distribution (|14p . It describes the distribution 
of thermodynamic fields on shortest time scales (but still 
^ to), over which the system can only explore a single 
metastable state for lack of time to cross a barrier to 
another state. 



E. Selfconsistency 

In order to determine A(x) self-consistently, we need 
to calculate the two point correlation function 



Q{x) = {SaSb)Ho,x^t=x 



dyP{x,y)m'^{x,y), (53) 



which together with Eq. ([7]) closes the self-consistency 
loop. 



Taking the derivative of Eq. (|53|) with respect to x 
and using the flow equations (|49l5ip . we obtain the more 
convenient form 



dQjx) 
dx 



dA{x) r , . (drn(x^ 

dyP{x,y)[-^ 



dx 



(54) 



The self-consistency equations can be solved itera- 
tively. At the n'^ stage, the iteration loop consists in 
computing 

Q„ = Q„({A„}) (55) 

from (j53p . and a new coupling matrix via 

/3A„+i = xr'(/3Q«) - l//3Qn- (56) 

The latter can be done conveniently in replica Fourier 
space, see App. |El For some details of the numerical 
implementation we refer the reader to App. |F] 



F. Marginal stability of the glass phase 

In this section, we show that throughout the low tem- 
perature phase the system is in a marginal state, pro- 
vided that the replica symmetry is continuously broken, 
or more precisely, provided that dA/dx is continuous and 
bounded as x — > 1. In other words, permanent marginal 
stability is guaranteed by the absence of step-like dis- 
continuities in A(a;) on the shortest time scales. This 
assumption is confirmed by the full numerical solution 
discussed in the following sectioni^ 

Let us rewrite the self-consistency equation ([7]) using 
the Fourier transform identity (|E15|) as 



f3{Qc-mx)) = Xi[rl{x)] 



where 



[Q]{x) = xQ{x) 



Qc ^ Q 



Qiy)dy, 



Q{x)dx. 



(57) 

(58) 
(59) 



We have also used the shorthand notation 

rl{x) = (3\A, - [A]{x)) + 1/(Q, - [Q]ix)), (60) 

recalling that rsc(a;) takes the place of an (average) 
Thomas-Fermi screening length on the time scale x. (In 
technical terms, r'^^ is the replica Fourier transform of 
the inverse polarizability operator —11"^, cf. (O.) How- 
ever, the term "screening length" should not be taken too 
literally, since we have seen that homogeneous screening 
breaks down below Tc- As we will see below r^c decreases 
with increasing time (decreasing x), reflecting that on 
long time scales the electron glass screens better on av- 
erage. 



11 



Taking the derivative of Eq. (|57p with respect to x, we 
find after a httle algebra 



1 



dQ{x)/dx 

dA/dx ~ ij[rl{x)]' 



Combining with Eq. (|54p we find the relation 



dyP{x,y) 



dm{x, y) 



dy 



(61) 



(62) 



which expresses the marginality of the full RSB solution 
at all scales x where A(a;) is continuous. 

Evaluating it for a; = 1, we find in particular 



1 



dy 



cosh"* (/??/) 



(63) 



extending the relation p4p from Tc to the low tempera- 
ture regime. 

As mentioned above, Eq. (|63p is equivalent to the con- 
dition for marginal stability on shortest time scales, or 
in technical terms, to the marginal stability of the free 
energy ([55]) with respect to replicon fluctuations^. Note 
that the above reasoning shows that any model with con- 
tinuous replica symmetry breaking is marginally stable 
throughout the low temperature phase, as was derived in 
a different way in Ref. 'Sg, too. 

Besides Eq. ([M]). an additional relation between the 
screening radius rsc(l) and the field distribution P{l,y) 
follows from the self-consistency of the spin-spin correla- 
tor, i.e., the compressibihty Eq. ([ST]) . 



dn 

dfi 



= /3(1-Q(1)) = 



pPil,y) 
TYTTTT- (64) 



cosh (/3y) 



Anticipating a large screening radius rsc(l) at low tem- 
perature, we can use the asymptotic expansion of xiy 
Eq. (O, and ip, Eq. to conclude that 



dy 



dy 



cosh^(/3y) 

cosh"^ (Py) 



1 



:(1)' 



/7rrsc(l) 



(65) 
(66) 



A consistent and regular low temperature solution of 
both equations requires a screening radius rsc(l) ^ l/T 
and a scaling form of the frozen field distribution in 
metastable states 



Pil,y)=T^^My/T), 



(67) 



with vE'3D(e) ^ ^ for e ^ 1. This will indeed be con- 
firmed by the full solution in the next section. 

Note that this scaling implies P(y) ^ y^ for y ^ T, 
which is exactly Efros and Shklovskii's prediction for the 
Coulomb gap in 3D. Simple arguments for Thomas Fermi 



screening in presence of this pseudogap also suggest an 
average screening on the scale rgc ^ i/T where thermal 
fiuctuations and Coulomb interactions compete. 

Interestingly, the saturation of the pseudogap expo- 
nent 2 = D — I comes out as a necessary consequence of 
the marginal stability condition, without requiring any 
further assumption. We may interpret this result as fol- 
lows. The fact that the dominant metastable states of 
the Coulomb glass are just marginally stable means that 
the system assumes the highest local density of states 
which is still just compatible with stability requirements. 
The marginal states exhaust all these degrees of freedom, 
which leads to the saturation of the Efros-Shklovskii ex- 
ponent. 

This result is interesting not only because it gives a jus- 
tification of the saturation, but also because the present 
theory takes into account the collective response of many 
neighbors to the presence of a given electron, and thus 
goes beyond a stability analysis with respect to single 
particle hops. Our result shows that such many particle 
effects do in fact not alter the gap exponent. However, 
we will see that it suppresses the coefficient of the power 
law describing the pseudogap in the density of states. 

Note that the exponent being 2 is not merely a co- 
incidence. In fact, in arbitrary dimensions D > 2 with 
J(r) ~ and J(fc) ~ l/k^~'^, an analysis along the 
same lines leads to P{l,y) = T^/"-^-^ oAy/T), with 
'^D,v{f) ^ e-^/'^^i for e 1, again in agreement with 
the stability arguments by Efros and Shklovskii. In the 
two-dimensional case with 1/r interactions, there are log- 
arithmic corrections (see the remark after Eq. p6)) ). 

The physical implications of the marginality of the 
glass will be discussed in Sec. |X] 



VI. LOW TEMPERATURE ANALYSIS 

Let us now analyze the low temperature behavior of 
the mean field equations in more detail. 

In order to solve the Parisi equations in the low tem- 
perature limit we introduce the rescaled variable b = x/T 
which appears in the anomalous response (I45p as the in- 
verse of an "effective temperature" T/x, replacing the 
real temperature for equilibrium response. Indeed, in 
Sec. IVIirB] we will review that the relation (|l5)) is actu- 
ally just the generalized fiuctuation-dissipation relation 
in disguise which is usually taken as a definition of "ef- 
fective temperature" . 

We anticipate the most important variations of Q{x) 
and A(x) on the scale of x ^ T, as is usually the case in 
mean field glasses^". Such a scaling has been observed 
in the recent numerical solution of the SK model at low 
temperature^, and has been exploited later to resolve 
the asymptotic T limit^^i^. 

The analysis in the previous section suggests that at 
low T an Efros-Shklovskii Coulomb gap develops, P(x = 
1, y) ~ i/^, and accordingly we expect typical overlaps 1 — 
Q ^ and screening radii r^c ^ 1/T. We will see below 
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that this scaHng not only holds on shortest time scales 
X = 1, but extends to larger time scales, too, provided the 
temperature is replaced by the "effective temperature" 
T/x = 1/6. In particular, we will find the overlap 1 — 
Q{x) and the average screening radius rsc{x; T) ^ 

b to decrease with increasing effective temperature, or 
equivalently, with increasing time. 



A. Regular flow equations 

The effective temperature sets a characteristic energy 
scale at time scale x. This suggests to rescale the local 
fields y by T/x, as was done in Ref. [sl for the SK model. 
We thus change variables to z = f{b)y = (c + b)y, where 
an arbitrary constant c = 0(1) is included to make the 
variable change regular for b — (3x 0. (When plot- 
ting results for p{z, b ~ (3) later, we will always revert 
to c = 0, however). Anticipating an Efros-Shklovskii- 



from the relations 



pseudogap P{y) " 
field distribution 



y , we further introduce the rescaled 



p(6,0,T) = b^P{x^bT,y^z/f{b)), (68) 

for which we expect p{z) ^ at large b = f3x and z ^ 1. 

In terms of the variables b and z, the magnetization 
and the rescaled field distribution obey the flow equations 



fib) 



dm 
'db 



dm 



dz 

bp{b) dK (dhn 
2~ d6 

-zf'ib)^+2p 
oz 

bp(b) dA fd^p 



(69) 



b d[m^] 
fib) dz 



(70) 



d[pm] \ 



db \dz^ f{b) dz J ' 



with boundary conditions that remain regular in the limit 
0, 



m{b:=:: /3,z) = tauh [z/ (/3) //?] , 
Pib^0,z) . exp{-[.//(0)]V2[I4^^ 



A(0)]} 



/2(0)V2^[V72 + A(0)] 



(71) 
.(72) 



The virtue of the rescalings is that /3 has disappeared 
from the flow equations, and only remains present in the 
boundary conditions for m. In the new variables the low 
temperature limit is completely regular, and focuses on 
the features of the physically active degrees of freedom 
described by z = 0(1). This will be exploited in the next 
section. 

The local spin-spin correlator (overlap) is obtained 



1 - Q(6 = P) 



1 



b'fib) y_c 

1 



dzp{b, z)[l — m{b, z)^] 
p{b = /3, z) 



P^fiP) J- 
dQ{b) _ 1 dA{b) 



dz 



cosh (z) 



db 



f(b) db 



dzp[b, z) 



dm{b, z) 



1 2 



dz 



b=f} 
(73) 

(74) 



From the first equation we expect 1 — Q{b) ~ b~^ . If 
we further assume little dependence on /3 at a fixed b 
we expect dQ/db ^ 5^^, which together with the second 
relation suggests dk/db ^ b~^ . This scaling will indeed 
be confirmed below at low but finite temperature, and 
will be used for the strict T = limit in Sec. IVIII 

The self-consistency condition Eq. ([F7|) . and its first 
derivative with respect to 6 = x/T, read 



Xi[rl{b^P)] 
d[rlib)] 



db 



ViM\ = -b 



/?[i-g(6 = /3)], 

dQ{b) 



db 



(75) 
(76) 



which allows us to obtain the average screening radius 
rsc(6) from the overlap function (3(6). We will show in 
Sec. IVIII Bl that Eq. ([75]) simply expresses Sompolinsky's 
anomalous response relation ([45]) at scale x = bT. 

Finally, we re-compute the coupling matrix A (6) by 
integrating the equations 



A(6 = 0) = Q(6 = 0)^[r2^(0)], 
dA(6) dQ{b) 



db 



db 



Mrtcib)], 



(77) 
(78) 



which closes the self-consistency loop A(6) 
rlib) ^ A(6). 



Qib) 



B. Results for the glass phase 

We have solved the self-consistency loop (|69ll78p nu- 
merically for temperatures down to T = 0.01, both for 
moderate and relatively strong disorder, W = 2, 10. De- 
tails of the implementation are given in App. |F1 

Our results confirm that the glass transition occurs 
in a continuous fashion at Tc, excluding a discontinuous 
breaking of replica symmetry, as is consistent with a sta- 
bility analysis in a Landau expansion around Tj,. The 
coupling function A{x) was found to remain regular ev- 
erywhere, ensuring marginality at all scales. 



1. Opening of the pseudogap 

In Figs. 131 E] the temperature evolution of the pseu- 
dogap below the glass transition is shown. Apart from 
the overall width of the field distributions, moderate and 
strong disorder {W ~ 2, 10) do not differ greatly. No- 
tice the fast opening of the pseudogap in the distribution 



FIG. 4; (Color online) Evolution of the pseudogap in the dis- 
tribution of the thermodynamic fields P{y) (top panel) and 
the instantaneous fields P{h) (bottom panel) for moderate 
disorder W — 2. The curves are plotted for the inverse tem- 
peratures P = 7.2, 7.6, 8, 10, 14, 20, 40, 60, 80, 100 (from top to 
bottom) below the glass transition {Pc = 7.149). Note the 
rapid opening of the pseudogap in P(y) just below the tran- 
sition. Most of the suppression of the instantaneous density 
of states is already formed as a plasma dip at the transition. 



FIG. 5: (Color online) Evolution of the gap in the distri- 
bution of the thermodynamic fields P{y) (top panel) and 
the instantaneous fields P{h) (bottom panel) in strong dis- 
order W = 10. The curves correspond to inverse temper- 
atures /? = 8.2, 8.4, 8.7, 9, 10, 14, 20, 40, 60, 80, 100 below the 
glass transition (/3c ~ 8.134). The salient features are similar 
as for moderate disorder, cf.. Fig. |4l 



of the thermodynamic fields. This translates into a fast 
drop of the compressibility below Tc- It may also affect 
transport properties to the extent that P{y) contains in- 
formation about the multiparticle excitations (electronic 
polarons) which are believed to carry hopping current^. 

On the other hand, the distribution of instantaneous 
fields displays a prominent plasma dip already at temper- 
atures substantially above Tc, as shown in Fig. [1] The 
further suppression of this dip affects mostly very low 
energies. 

Both for P{h) and P{y) only that low energy regime 
is expected to exhibit a universal ES pseudogap. For 
better comparison of the temperature evolution of the 
two distributions, we display them together in Fig. [6] for 

T > Tc, T w and T < T^. 



2. Compressibility 

In the main panel of Fig. [7] we plot the compressibility 
(j64p as function of T. The higher order non-analyticity 
at the transition temperature Tc is hard to discern. It 
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FIG. 6: (Color online) The distribution of instantaneous fields 
h (solid lines) and thermodynamic fields y (dashed lines) for 
temperatures well above (/3 = 2, red), close to (/3 = 7.2, green) 
and well below (/3 = 100, blue) the glass transition at mod- 
erate disorder W = 2. Notice the presence of a pronounced 
plasma dip in P{h) even at high temperature. 
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would in principle show up more markedly in a plot of 
the non-linear capacitance, which is the analog of the 
non-linear susceptibility of spin glasses. 

Similarly as in spin glasses, the low temperature, zero- 
field cooled compressibility is strongly suppressed as 

Kc^nW^ ^m-Qa)]-T\ (79) 

reflecting the lack of short-time ergodicity of the system. 
On the other hand, the equilibrium, or field cooled, pre- 
diction for the compressibility is 

pi 

dx [1 - Q{x)] const. (80) 

Jo 

This observable does not feel the Coulomb gap and is 
thus much less suppressed. The field-cooled compress- 
ibility is only accessible by raising the temperature back 
above Tc and quenching the system in the new field, or 
alternatively, in the ultra-long time limit at a constant 
low temperature. The unstable, replica-symmetric high 
temperature solution similarly predicts a constant sus- 
ceptibility as T ^ 0. 

3. Tunneling 

Let us briefly discuss tunneling from a broad junction 
into a strongly insulating classical electron glasses, as 
studied experimentally in Ref. [3^ . The broad junction 
assures self- averaging of the tunneling current. In the 
insulating limit, and for voltages well below the temper- 
ature, eV < ksT, the tunneling current is dominated by 
the hopping of electrons from the probe to empty sites 
in the electron glass. Such hops are usually not accom- 
panied by any further rearrangements of other electrons, 
since those are exponentially suppressed in terms of the 
action. 

In contrast to the compressibility, tunneling is thus 
controlled by the distribution of instantaneous fields, 
P{h), pT|) . In the linear response regime the differen- 
tial tunneling conductance is proportional to 

G,r,.^,K^O)c</..-i^, ,81, 

which is plotted in the inset of Figs. [71 

C. Scaling and universality 

The low temperature solution confirms the anticipated 
scalings dQ/db ^ b~^, dA/db ^ b^^ and rgc ^ b for large 
f'l ^ b < (3. In this regime these functions, at fixed b, 
are nearly independent of temperature and disorder, as 
illustrated in Figs. [5] and [51 As a result, physical quan- 
tities such as the compressibility, or the density of states 
of thermally active degrees of freedom turn out to be re- 
markably universal, that is, independent of disorder and 
obeying simple scaling laws as a function of temperature. 




T 



FIG. 7: (Color online) Main panel: Zero-field cooled com- 
pressibility Kc (ZFC, blue line), equilibrium field cooled com- 
pressibility (FC, green line) and the replica symmetric result 
(RS, red line) for W = 2. The extrapolation of the RS pre- 
diction below Tc is unphysical, but closely follows Kq'' , as one 
may expect. Inset: Tunneling conductance G as predicted 
by Eq. (|8ip in the glass phase (RSB). The unphysical replica 
symmetric theory (RS) would predict a strong suppression of 
G due to a constant Onsager term and an ensuing hard gap 
in Prs(/i). 
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FIG. 8: (Color online) The coupling function A(6) for /3 = 
20, 40, 60, 80, 100, and disorder W = 2. The result is nearly 
independent of /3 as long as (3c <^ b < 0.7/3. 



Furthermore, we find the scaled field distribution func- 
tion p{b = /3, z) to approach an essentially temperature 
independent form at low T (see the top panel of Fig. [T5)) . 
This confirms our prediction (|67p for the scaling form of 

The reason for the asymptotic universality at low T 
lies in the appearance of an attractive fixed point of the 
fiow equations, as discussed in the next section. 
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consistency equations (I75II78P to 




FIG. 9: (Color online) The coupling function A(&) at low tem- 
perature T = 0.01, for different disorders, W — 2, 10. In the 
regime b S> /9c, the coupling function is universal, i.e., essen- 
tially independent of the bare disorder. 



D. Break point and plateau 

On short timescales {x ^ 1, b (3), the sole depen- 
dence on x/T is violated. The derivative of the cou- 
phng function, b^dA/db, drops to rather sharply around 
xp{T) = bpiT)/(3 < 1, implying that both A(a;) and 
Q{x) are constant beyond xp(T). This plateau in the 
coupling and overlap function appears below T^, with 
xp{T) decreasing to a saturation value xp{T ^ 0) ~ 0.78 
at sufficiently low temperature. The low temperature 
limit of this feature and its physical implications will be 
analyzed in the following two sections. 



VII. LOW TEMPERATURE FIXED POINT 

A. Flow equations for T = Q 

Let us now analyze the flow equations in the regime 
of very low effective temperatures, T/x — l/b ^ Tc- In 
particular we will be interested in the zero temperature 
limit at finite x where x 3> T/Tc is satisfied. For this 
analysis we can use the asymptotic low energy scalings 

of Eqs. mm 



2-2/q 



(82) 



where we allowed for a more general asymptotics for ip so 
as to cover arbitrary space dimensions D, with power law 
repulsions 1 /r''. Note that the parameters A and a derive 
from the low energy asymptotics of the interactions Jy , 
and are thus independent of lattice details. For Coulomb 
glasses (i^ = 1), a is the Efros-Shklovskii exponent a = 
D ~ 1, and in particular A^n = ^/^^, cf., Eq. (fT8| . 

The above asymptotic forms simplify the self- 



A(x) = A[r,c(x)]2-2/"Q(x), 
[riKx) 



= l3xQ{x), 



(83) 
(84) 



where dots denote derivatives with respect to x. 

The low temperature analysis of the SK model^^ leads 
to completely analogous equations, the only trace of its 
specific interactions being the values A ^ 1, a = 1 
(since V'SK = !)■ The low energy asymptotics of the SK 
model is thus very similar to a 2D electron glass, apart 
from logarithmic corrections in the latter which lead to 
V'2D[r2,] «27rln[rsc(a;)2]. 

In order to obtain a regular low temperature limit of 
the flow equations, we rescale the physical quantities sim- 
ilarly as in the previous section (without changing the 
variables x — > 6, however). 



z 

m{x,y) 
P{x,y) 
A(x) 
Qix) 
lAsc(a^) 



(3xy, 
m{x, z), 
(/3a;)"Xx,z), 
5l3{(3x)-^--^\{x), 

CP{pxy^'\{x), 

iPxy^ix). 



(85) 
(86) 
(87) 



(89) 
(90) 



The exponents 77, 5, P are chosen in such a way that any 
explicit singular dependence on (5 is eliminated from the 
equations. The form of the flow equations for ni and P 
and Eqs. (|83l84p together require 



C-1 



5 = 2, (91) 
D — 1 Electron glass, 



1 



SKmodel. 



(92) 



With these notations the flow equations become: 



xm — — A(x) (to" -|- 2mm!) — zm' , (93) 
xp = X{x) {p" - 2{pm)') - zp' + ap, (94) 

where primes denote partial derivatives d/dz. The self- 
consistency equations (1831184^ take the form 

\{x) = |(l + a)g(x)e(i/"-i)(x), (95) 



xS^{x) — aS,{x) — {I + a)q{x) , 
while Eqs. (|73l74p translate into 



Cq{x) = 6X{x) I dzpm'^, 



= 1) 



dz ■ 



cosh (z) 



(96) 

(97) 
(98) 



Notice that for the SK model {a — 1) Eq. expresses 
simply the identity A(a;) = q{x), while Eq. (|96p is not 
needed to obtain a self-consistent solution. 



B. Fixed point 



The system of Eqs. (|93l94p looks like a set of func- 
tional renormalization group (FRG) equations as they 
appear in studies of vortex lattices, magnetic interfaces 
and wetting problems in the presence of randomnes a^'^'^^ . 
In the present case, the "time variable" x takes the role 
of the renormalization scale, m and p being analogous to 
functional running coupling constants. However, a for- 
mal difference with standard FRG equations consists in 
the explicit x-dependence of the " beta-function" through 
X{x) which is itself determined by the flow. If X{x) 
were constant, the flow equations would constitute an 
autonomous system of partial differential equations (i.e., 
with a right hand side independent of x). Such systems 
commonly possess fixed points, such as many standard 
renormalization group flows. 

It turns out that despite the remaining x-dependence 
the flow is attracted to an asymptotic fixed point in the 
intermediate regime 1 3> a; ^ T/Tc, which controls uni- 
versal low energy physics in a similar manner as in stan- 
dard renormalization group flows. The flow of the mag- 
netization and field distribution functions to respective 
fixed points is illustrated in Figs. [TOl and [TT] 

The attractiveness of the fixed point wipes out the 
trace of the bare disorder W entering the flow equations 
via the boundary conditions at cc = 0, that is, at high ef- 
fective temperature and energy scales. The flow to lower 
effective temperatures renormalizes the field distribution 
p{z) to a universal form, at least for fields describing ac- 
tive degrees of freedom, y ^ Tctf or z — 0{1)). On the 
other hand, the high energy tails in the unsealed distri- 
bution P{y) basically retain the structure set by the bare 
disorder. 

Despite the above suggestive similarities with genuine 
functional renormalization group flows, we do not have 
- at present - a thorough understanding of this analogy, 
even though we believe that there is a deeper connection 
with FRG approaches. 

On the putative fixed point, the derivatives x-^ van- 
ish. The fixed point functions m*,p* depend only on the 
variable z, satisfying the equations 



A* (m*" + 2m*m*') + zm*' ^ 0, 

A* (p*" - 2{p*m*y) - zp*' + ap* = 0, 

(q* = SX* f dzp*{m*'f, 



aC = {l + a)q\ 



with the self-consistent value for A* 
A* 



2 ^ ^ 



(99) 
(100) 

(101) 

(102) 



(103) 



Aa f 1 + a . 



I/O 




FIG. 10: (Color online) Illustration of the temporal flow 
of m{x,z) for increasing Sompolinsky times x = b/P — 
1,0.9, 0.8, ...0.2 for (5 = 100 and W = 2. An intermediate 
fixed point is clearly approached near x — 0.5. The asymp- 
totic prediction for the fixed point is plotted as a dashed line 
(black), for which we used an artificial power law fitting ipso 
beyond its universal low energy regime. 




FIG. 11; (Color online) Illustration of the temporal fiow 
of p{x,z) for increasing Sompolinsky times x — b/ (3 — 
1.0,0.9, . . .0.2 for /3 = 100 and W = 2. The apparent fixed 
point occurring around x = 0.5 comes close to the fixed point 
predicted with an artificial power law which fits beyond 
its universal low energy regime (dashed black line). The re- 
maining small discrepancy is due to the approximate nature 
of the power law fit. 



conditions 



TO*(0) =p*'(0) = 0, 

m*{z 00) —>■ sign(z). 



p*{z — !■ 00) cx 




The fixed point functions are subject to the boundary 



(104) 
(105) 

, (106) 



where Ha is the Hermite polynomial of order a. The last 
condition follows from the solution of the linear Eq. (|100[) 
at large z where m* — sign(z). This leaves open a con- 
stant of proportionality which has to be fixed eventually 
by the self-consistency constraint (|103p . 
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Eqs. (|99ll00p constitute a self-consistency problem for 
A*. They have a unique solution which can easily be 
found numerically. For the 3D Coulomb glass we find the 
fixed point values A* = 1.1230287..., q* = 0.26763357..., 
^* = 0.40145035.... For comparison, in the SK model one 
finds \* = q* = .4108021... and ^ = 0.8216042.... 

C. Stability 

The analysis and even the meaning of the attractive- 
ness of the fixed point is not as straightforward as in 
standard renormalization flows, because of the non-local 
coupling due to the dependence of A(a;) on the flow itself. 
However, we can analyze deviations from the fixed point 
of the form 

m{x,z) — m*(z) + mi(z)a;'', (107) 
p{x,z) = p*(z)+pi(z)a:^ (108) 
A(x) = A^ + Aix"^, (109) 

such that the consistency equations (|95ll98p are still 
obeyed. The functions mi, pi and Ai will in general de- 
pend on the exponent i/. We found that both for the 3D 
Coulomb glass and the SK model such linearized devia- 
tions require a positive v > fmin > 0. This suggests that 
the fixed point is attractive for decreasing x. Indeed, the 
full solution approaches the fixed point more and more 
closely at small x 3> T/Tc, as is clearly observed when 
one solves the T = flow equations numerically. How- 
ever, it is not entirely clear to us what physical meaning 
should be ascribed to this kind of attractiveness. 

Once the existence of the fixed point and its attraction 
for decreasing x has been established, we can take p*{z) 
as an initial condition for very small x and forward inte- 
grate towards x — 1. This last stage may be viewed as 
a renormalization flow of the fleld distribution to its low 
energy effective form, which attains an invariant form on 
intermediate time scales. 



D. Low T solution 

At sufficiently low temperature, the fixed point con- 
trols the fiow of p{x, z) for z — 0{1) for all x ^ T, in the 
sense that one can take p*{z) as an initial condition at 
a; » T. In the fimit T ^ 0, one then solves Eqs. (|93II98|) 
for all X in an iterative manner analogous to the proce- 
dure described in the previous section. 

At this stage, one may worry about the influence of 
boundary conditions for p{x, z) at large z. However, it 
turns out that the fixed point is very stable and strongly 
attractive in the considered models, such that the precise 
form of the boundary conditions has very little influence 
on the resulting flow. Their effect seems to be exponen- 
tially small in the value Zhc where the boundary condi- 
tions are imposed. A rigorous way to ensure small errors 
from boundary conditions at finite but low temperature 



is described in App. [F] The latter allowed us to assess 
the good quality of the results we obtained with simple 
boundary conditions resulting from forward integrating 
the fiow for p in the regime z 3> 1 where m — sign(z), 

p{x, z) = z" + 772"-i(2A* - xtJi{x)) + ©(z"-^). (110) 

Here, Ah{x) is given by 

Ahix) = 2 r dx' ^^'^'\^^* . (Ill) 
Jo ^ 



E. Scaling of the density of states 

The fixed point distribution p*{z) governs the low tem- 
perature behavior of the thermodynamic field distribu- 
tion of metastable states, p{x — 1, z) — P{x — l,y — 
zT)/T'^ in the sense that it leads to a universal scaling 
function p(a: = 1, z) ^ \I'3d(z), which is shown in Fig.fT2l 
While the fiow from the fixed point to cc = 1 (relating p* 
and ^3d) modifies the precise functional shape on the 
scale of z '-^ 1, it preserves the large field asymptotics, 
see Fig. 111! In particular, we may extract the asymptotic 
value for the prefactor of the parabolic pseudogap from 
the fixed point function at large z where to very good ap- 
proximation we have p(x = 1, z) ~ p*(z) = Cz^ -I- 0(z), 
with C = 0.02604. This implies P{y) = Cy^ + 0{y) for 
T ^ y ^ Tc- Translating the local field to the energy to 
introduce additional particles (or to flip the local spin), 
E = 2y, and restoring the lattice units e'^/ini and £ for 
energy and length, respectively, we flnd for the single 
particle density of states 

p{E) « SC'^E^ = 0.2083^E'^, T <^ E 1^.(112) 

e° e" 

The coefficient of the parabolic term is significantly 
smaller than the value 3/tt « 0.955 predicted on the basis 
of the minimal stability requirement of a low energy con- 
figuration with respect to single particle hop o^^'^^ . This 
can be ascribed to the fact that our mean field theory 
takes into account the back reaction of many particles, 
which should indeed lead to a more stringent restriction 
for the local density of states. It is interesting to note 
that an estimateSS. based on a simple function P{y) and 
imposing marginality, leads to 8C « 0.204 very close to 
the exact value in Eq. ^T^^. 



F. Corrections to scaling 

So far we have concentrated on the extreme low en- 
ergy limit where we approximated i/'('"sc) by the univer- 
sal term in Eq. p8)) . In order for this term to domi- 
nate we need r^^, :s> (2Ciatt)^/7i'. From the knowledge of 
the fixed point asymptotics r^^ « {/3x)^/^* for x < xp, 
we estimate that we need to go to temperatures below 
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FIG. 12: Universal, asymptotic low temperature scaling func- 
tion *3d(2) = P{l,y = zT)/T^ = p(l,«) (full line) and the 
fixed point function p*{z) (dashed line). Both have the same 
curvature at large z. 



T <^T* = v'tt/C :rp/2|Ciatt| « 0.12 = Qme^/ni to en- 
sure universal asymptotic fixed point behavior. At tem- 
peratures and energies of the order of or larger than T*, 
the non-universal, lattice dependent corrections (|10p in- 
fluence the physical observables, in particular the density 
of states. 

The effect of those corrections are clearly visible in the 
rather slow approach of our finite temperature data for 
p{z,x = 1) to a universal function \E'3d(z), see the top 
panel of Fig.[T31 This can be traced back to the deviation 
of V^SD from its low energy asymptotics for rgc 2> 1. The 
convergence of the finite temperature data looks much 
more convincing if we repeat the above analysis, but ap- 
proximate '03d(?'sc) with an artificial power law Ar^ 
which fits the function best in the relevant regime of pa- 
rameters Tsc. For the temperature range 30 < /3 < 100, 
an exponent a ~ 2.34 turns out to provide a reasonably 
good fit to i/jsD , and to produce a much better scaling of 
p{z) = /3"P(y = 2;//3), as shown in the bottom panel of 
Fig. [T31 As a consequence, it looks as if the density of 
states were governed by a larger exponent a at energies 
E > T*. 

The apparent fixed point in the fiow of average mag- 
netizations and field distributions as a function of Som- 
polinsky time x is well captured by the fixed point theory 
applied to the above artificial power law, as illustrated in 
Figs. [TO] and [TT] However, the genuinely universal low 
temperature fixed point predicted by Eqs. (|99ll00p only 
applies to the very low energy and temperature regime. 

A similar effective power law of the pseudogap was 
also found in numerical simulations of 3D Coulomb sys- 
tems. The energies for which the density of states can be 
sampled reliably lies systematically above T* , because 
of finite size limitations. For this non- universal regime 
Mobius et ali^ reported exponents a of the order of 
2.2—2.4 at the lowest reliable energies, in good agreement 
with our mean field prediction. 




FIG. 13: (Color online) Scaled distribution of thermodynamic 
fields as a function of temperature. The curves correspond 
to = 2 and /3 = 20,40,60,80,100 (from top to bottom). 
Top panel: Scaling with the Efros-Shklovskii-exponent a = 2, 
p(l,z) = P{l,y = zT)/T^. The inset demonstrates that the 
differences AP = |Pvy=2 — Pw=io| due to the strength of 
disorder tend to zero in the low T limit, as expected from 
the universality induced by the fixed point. Bottom panel: 
Improved scaling with an artificial power a = 2.34, p(l, z) = 

p(i,y = 2r)/r". 



G. The plateau and the overlap distribution at 

r = 

The T ^ limit of the self-consistent coupling func- 
tion X{x) is shown in Fig. 1141 Starting at its fixed point 
value A* at small x, it slightly decreases with increas- 
ing X, and finally vanishes rather abruptly at the plateau 
threshold xp ^ xp{T ^ 0) ^ 0.78. The function q{x) 
obviously exhibits an analogous behavior. 

For X > xp, X{x) = q{x) = which implies that 
the overlap function Q{x) has a plateau for xp < x < 
1. This has interesting consequences: According to 
Parisi's interpretation of the overlap function, P{Q) — 
[dQ / dx]"^] I , describes the distribution of overlaps 

among metastable states, if sampled according to their 
Boltzmann weight. The plateau implies that the low- 
est lying state (the ground state valley) takes a finite 
weight of the total Boltzmann weight, contributing a 5- 
peak at the Edwards- Anderson overlap QeA: correspond- 
ing to the self-overlap of typical low-lying states. Higher- 
lying states are separated from the lowest state at least 



19 




FIG. 14: T = solution of the coupling function X{x) — 
{Px)^A{x)/2f3 for the 3D electron glass. The fixed point 
controls the value X{x 0) = A*. The drop to around 
xp « 0.78 indicates a plateau in the coupling and overlap 
function. 



0.75 



<T3 




0.25 



FIG. 15; (Color online) Equilibrium distribution of the 
rescaled overlap distance, d = {1 — Q) /T^ , in the 3D electron 
glass. There is a peak of strength 1 — xp at dsA = 0.8609 
corresponding to the intrastate (Edwards- Anderson) overlap, 
and an asymptotic power law tail P{d) ~ (|-(2+a)/(i+") ^ 



VIII. PHYSICAL INTERPRETATION OF THE 
ASYMPTOTIC FIXED POINT 

Apart from providing accurate numerical values of low 
temperature characteristics, the asymptotic fixed point 
in the flow equations has several interesting physical con- 
sequences. 



A. Dynamical self-similarity 

Consider the distribution of local magnetizations aver- 
aged over a time scale tx, p{m; x). On intermediate time 
scales X, with T/Tc <^ x < xp, the fixed point attracts 
the flow and we flnd to a good approximation 



p(to; x) 



dy P{x, y)5{m — •m{x, y)) 
1 



{Pxf 
1 



dzp*{z)S{m — m*{z)) 



{f3xf 

with the invariant form factor 



p^m)^T,%{x)p*iml (114) 



p*{m) 



p*iz) 



m*'{z) 



(115) 



m—m* (z) 



This applies for all |1 — m| such that z — [m*] ^ (to) is 
not too large (|1 — to| not exponentially small), such that 
the flxed point function indeed applies. 

Notice that the dynamics merely increases the num- 
ber of active degrees of freedom (by the third power of 
the increasing factor Tcs/T , as seems natural in presence 
of a parabolic pseudogap), without changing the overall 
shape of the distribution except in the tails. This result 
indicates an interesting self-similarity of the dynamics. 
We may actually speculate that after a suitable coarse- 
graining of time the dynamics will look essentially the 
same on all time scales, except that the actual number 
of active degrees of freedom gradually increases. 



by 0{T) in energy. Similar features^ are known for the 
SK model2^^. 



The scaling ([M)) suggests to introduce the reduced 
phase space distance d = (1 — Q)/T^, with the intra- 
state (Edwards-Anderson) distance d^A = 0.8609. The 
asymptotic fixed point q{x) — *■ q* implies a power law 
decay of the tails of the overlap distribution at low T, 



P{d) = {l-xp)S{d-dEA) + xpPtM, (113) 



where Pts^nid) ^ d (2+a)/(i+a)_ xhe low temperature 
limit of this distribution is plotted in Fig. [T51 



B. Generalized fluctuation-dissipation relation 

1. Exact relation for global properties 

In equilibrium, the classical fluctuation-dissipation 
theorem (FDT) relates the integrated response function 
x{t,t') and the correlation function C{t,t') via 



(116) 



In equilibrium, time translation invariance holds and 
both functions only depend on t — t' . However, in a 
glass, time translation invariance is broken because the 
system remains out-of-equilibrium and "ages". That is, 
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over longer and longer times it explores growing portions 
of phase space. In this case, the relation (|116p no longer 
holds on large time scales, but is replaced by 



gration to show 



R{t,t') = (3x{C) 



dC{t,t') 
dt' 



(117) 



The extra factor x{C) < 1 only depends on the corre- 
lation C, provided both times are large, t,t' ^ to. The 
temperature in the FDT relation simply is replaced by 
Tcsit,t') = T/x{C{t,t')), which, at fixed initial time t' , 
may serve as a measure of the observation time t, as we 
used implicitly all along this paper. 

The above relation (|117[) actually defines an " effective 

temperature" as the ratio '^^glf ^ /R{t, t'). 

Out of equilibrium dynamical behavior such as de- 
scribed above appears explicitly in the mean field solu- 
tion of the Langevin dynamics of various spin glass mod- 
glg55,56^ For the case of glasses with continuous replica 
symmetry breaking, one can obtain the generalized FDT 
relation (|117p directly from Sompolinsky's anomalous re- 
sponse Ansatz ps)) . In fact, pi7p is just a simple rewrit- 
ing of the latter, if we note that C{x) appearing in the 
FDT relation is the same as Q{x) in Sompolinsky's dy- 
namics. Therefore, for such mean field models, the x in 
pi7p and the one appearing in the replica and Sompolin- 
sky formalisms are in fact identical. 

Even though these arguments may suffice to show that 
the generalization of the FDT relation is automatically 
built into Parisi's replica symmetry breaking scheme, it is 
useful to derive the relation directly from the flow equa- 
tions. This will illustrate the meaning of the quanti- 
ties m(x, y) and P{x, y) and the formalism introduced to 
solve the mean field equations. 

In the spirit of Sompolinsky's approach, the integrated 
response function x(*:^') = /(/ R{t,t")dt" takes the nat- 
ural expression 



Pix,yf-^dy 
dy 



(118) 



where x is associated to the pair of initial and observation 
times {t' ,t) in the sense explained above. In the replica 
formalism one anticipates that this response is given by 
f3{saSh)- By inspection of Eqs. (I6l7p . this suggests the 
relation 



(119) 



which we will confirm shortly. 

The time averaged correlation function is similarly 
given by 



C{t,t')= P{x,y)m\x,y)dy=Q{x). 



(120) 



Taking the time or x-derivative of both equations 
pi8ll20p , one easily establishes the generalized FDT re- 
lation pi7p by using the flow equations and partial inte- 



(121) 



ax 



To prove the second relation, we have used Eq. ([M]) . 

Further, a comparison with Eqs. ([57|) or (|76|) estab- 
lishes the equivalence conjectured in Eq. ()119|) . 



2. Local interpretation of T^s 

The above relation is an exact global relation which 
makes a statement about the correlations and response 
averaged over the whole system. On the other hand, 
the asymptotic fixed point allows us to establish a local, 
though only asymptotically exact, counterpart. 

Consider the average local magnetization of a site i 
seeing a frozen field yi over a time scale with the same 
restrictions on x as earlier. In that regime, it is governed 
by the fixed point 



(sj)t^ « m*[yi/Tcs[x)]. 



(122) 



Moreover, the time averaged response to an external field 
/loxt applied over a comparable time scale is governed by 
Toff instead of T 



1 



^{mi)t, __ 

d/lext Toff (2;) 



{m*)'[y/T,s{x)], (123) 



where m* substitutes for the familiar tanh as a response 
function. 



IX. AGING 

A. Physical meaning of the plateau 

An appealing static interpretation of Parisi's Ansatz 
was found by Mezard et ali^. They rederived the so- 
lution within a cavity approach, by assuming a hierar- 
chical organization of states within clusters within even 
larger valleys etc. An essential ingredient is the distribu- 
tion of the free energies of states relative to the cluster 
they belong to, the distribution of cluster energies rela- 
tive to their valley, and so on. In order to recover Parisi's 
RSB scheme for the SK model in a if step RSB approx- 
imation, one has to suppose exponential distributions 
p{F^^'') (X exp [xK+i-tPP'^^^], where the f are the free 
energies pertaining to a valley on the i*'' level of the hier- 
archy {i — 1, . . . , K). The parameters 1 > xk > ■ ■ ■ > xi 
extremizing the free energy turn out to be identical to 
the optimal size of replica blocks in a if-step approxima- 
tion. This gives a further interpretation of the variable 
X as describing the energy landscape at a certain level of 
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the hierarchy of valleys (which is explored by the system 
on the Sompolinsky time scale corresponding to x) . 

For experimental purposes, one would like to know how 
the immediate environment of a state is organized, since 
this is the local landscape being usually explored in glassy 
response measurements. According to the above and the 
ultrametric picture of phase space organization, the free 
energies of states close to a given metastable state are dis- 
tributed exponentially, and parametrized by the largest 
parameter xk appearing in the above hierarchy. This 
picture carries through to the case of continuous RSB 
where this largest value is given by the plateau threshold 
xp = lim^^oo xk- 

B. The trap model and simple aging 

The information about the local neighborhood of a 
metastable state can be combined with Bouchaud's trap 
model^i22 in order to make a non-trivial prediction about 
the aging behavior of the electron glass. 

The trap model assumes the dynamics of the glass to 
be a random trajectory among the set of states within 
one cluster of Parisi's hierarchy. The escape from a given 
state ("trap") j is assumed to be Poisson distributed with 
an average escape time tj = exp[PFj], where Fj is the 
free energy barrier associated with the j"^ trap. Tak- 
ing the above cavity interpretation of Parisi's ansatz as 
a motivation, the Fj are taken to be random variables 
with a probability density p{F) ~ exp[(3xFj], so far with 
a free parameter x. Finally, a real system's dynamics 
is expected to perform a natural ensemble average over 
such trap models, since different regions in space can be 
thought of as roughly independent trap models. 

For an exponential barrier distribution, the trap model 
dynamics depends sensitively on the parameter x. For 
X > 1 the system is ergodic, sampling all of traps in 
a time proportional to their total number. However, if 
X <1 aging effects occur, reflecting that the dynamics is 
weakly non-ergodic: While the system is not stuck in any 
given trap forever, it is nevertheless unable to explore the 
whole phase space in a finite time since the mean dwelling 
time in a trap diverges, (r) =00. 

To quantify these effects, one can look at a system 
that evolved from random initial conditions for a time 
tyj. The probability to find the system still in the same 
trap at time t as at tw has been calculated^^ to decay 
asymptotically like 

C{t;t^)r^[t^/tf, (124) 

where t and are measured from the same origin of 
time. 

The problem of connecting the phenomenological trap 
model with quantitative models for glasses lies usually in 
the lack of prediction for the aging exponent x. However, 
from our low temperature solution of the electron glass 
problem, we obtain such a prediction as x = a;p(T — s- 
0) ~ 0.78. Hereby, we think of the dynamics as exploring 
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FIG. 16; (Color online) Schematic representation of the en- 
ergy landscape assumed in the trap model. The depth of the 
traps is exponentially distributed, deeper traps being rarer 
than shallow ones. To make contact with the replica solution 
we have to assume that the states in a cluster have a common 
"ceiling" or threshold level beyond which the system easily 
escapes and becomes trapped in any other state within the 
cluster. 

only the lowest level of hierarchy in an ultrametric or- 
ganization of states. Moreover, we make the assumption 
that the states within a valley have a common threshold, 
as depicted in Fig. This allows us to identify the free 
energies F*^^-* of the cavity approach with the free energy 
barriers of an effective trap model, whose parameter x 
corresponds to the plateau value xp, the first (largest) x 
for which there is a non-trivial (non-constant) structure 
of the Parisi matrix Q{x). 

The above prediction may be relevant for aging experi- 
ments in electronic glasses. Indeed, recent measurements 
of electronic aging in 2D indium-oxide filmsi^ found con- 
ductance relaxations which obey "simple aging", i.e., a 
scaling of the relaxation with t/t^, with an asymptotic 
decay that can well be fitted to the power law (|124p with 
exponents x = 0.85 — 0.9. 

Similar experiments on granular aluminuroi^ seem to 
be consistent with an exponent a; = 1. A possible reason 
may be the different nature of the Coulomb interactions 
in such a system, which is suggested by the unusual trans- 
port behavior in the hopping regimei^ 



X. DISCUSSION 
A. Experimental implications of the glassy state 

The glass transition should in principle be observable 
as a divergence of the non-linear capacitance. In practice, 
this will prove rather difficult to measure, even though 
attempt sin this direction were madeS*-. It may be eas- 
ier to look for other indications of the glass phase, such 
as out-of-equilibrium behavior or manifestations of the 
glassy criticality predicted for the whole low tempera- 
ture phase. The criticality and marginal stability of typ- 
ical metastable states arises from their approximate de- 
generacy with many other nearby valleys in the energy 
landscape. The latter leads to a large number of collec- 
tive soft modes (low energy rearrangements) around local 
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minima and saddles in configuration space. 

To the extent that this prediction is not an artifact of 
the mean field approximation, it should have important 
physical consequences for the Coulomb glass. We ex- 
pect the marginality to be reflected in the original lattice 
problem as a permanent criticality of the low tempera- 
ture state, in the sense that the square of the connected 
two point correlator (|4ip retains a power law decay even 
below Tc- This implies in particular that screening is not 
exponential. Instead the charge correlations decay slowly 
and even have a random signs at far enough distances. 
Similar behavior was indeed observed in numerical sim- 
ulationsI2, at T = 0. 

We infer the above scenario from a closely analogous 
situation in the SK spin glass, where the local approxi- 
mation and Parisi's solution are known to be exact. In 
the same way as for the Coulomb glass, one can prove 
the marginal stability of the low temperature phase for 
all T < Tc in the single site model. On the other hand, 
it is known from the analysis of the Thouless- Anderson- 
Palmer equations for all local magnetizations^- m^, that 
the marginal stability of the single site approximation is 
directly related to an extensive number of gapless modes 
in the excitation spectrum around typical metastable 
minima^. We conjecture that a similar connection holds 
between the marginality of the local approximation and 
the criticality of the lattice Coulomb glass. If this is true, 
the presence of a large number of soft collective excita- 
tions will be a crucial ingredient to the understanding of 
glassy dynamics and relaxation. In the quantum version 
of the Coulomb glass, the marginality is closely related 
to the closure of the gap in the spectral function^**, and 
one expects the presence of the associated collective soft 
modes to have an important impact on conductivity and 
DC variable range hopping in particular. 



C. Limitation of the mean field approximation 

In the mean field approach we neglect certain corre- 
lations among the neighbors of a given "cavity" site. 
More precisely, the response of neighboring sites to the 
presence of a cavity spin is treated as essentially inde- 
pendent, the only correlations retained being ring dia- 
grams. This approximation does not include degrees of 
freedom such as electronic dipole a^^i^^ . These are pairs 
of sites, one occupied the other empty, with individual 
local fields I/11.2I 3> T, but with an excitation energy 
hi + h2 — Jij < T for swapping the particle hole pair. 
Such dipoles can be defined unambiguously only at tem- 
peratures sufficiently below T^. 

The interaction of single (active) sites with such de- 
grees of freedom is not, or not entirely contained in the 
single site approximation. However, it is possible that 
the behavior of single sites is actually dominated by the 
interaction with other active sites, rather than by the 
interaction with dipoles, which is assumed in our mean 
field approach. Indeed, we obtain very reasonable results 
for the single particle density of states. However, a clear 
indication for the neglect of dipolar degrees of freedom 
in the mean field theory is found in the specific heat. 
While the single site model would predict Cy ~ T^, it 
is known that dipoles dominate the specific heat of the 
full lattice model, with an almost linear temperature de- 
pendenceSii^iSa. 

It is a very difficult open problem to build a self- 
consistent mean field theory which would include such 
non-local objects. They must be expected to be dynam- 
ically generated as the system settles into low energy 
metastable states, their location and orientation depend- 
ing on the specific metastable state under consideration, 
which makes the problem very difficult to tackle. 



B. Extension to 2D electron glasses 

In two-dimensional electron glasses with l/r interac- 
tions, the mean field approximation is less well justified 
than in 3D due to the reduced number of effective neigh- 
bors. If we nevertheless carry out analogous approxi- 
mations as in 3D, we obtain a similar single site prob- 
lem, though with a function '!/'2d['"sc] ~ 27rln [rsc(a;)^]. 
Because of the logarithmic corrections, there is no gen- 
uine fixed point in the flow equations as in 3D. How- 
ever, in order to qualitatively describe the distribution of 
fields around a given T ^ y <C 1, we can approximate 
■02D ~ const. K, A K, 27rln(l/y), which corresponds to 
an exponent a = 1, reducing the single site model to an 
effective SK model with random fields. From the solu- 
tion of the latteriSi^ we know that = 0.30|z|/^, 
and we thus expect the 2D electron glass to exhibit a 
quasilinear pseudogap P{y) ~ 0.3|?/|/27rln (l/|y|) at low 
energies, with logarithmic corrections suppressing small 
fields. 



D. Comparison with numerical studies 

1. Glass transition 

So far numerical studies for the model ([T]) with on- 
site disorder have remained inconclusive as to the exis- 
tence of a glass transition, even though a transition was 
claimed to exist in a 3D system with positional disor- 
der but no on-site randomness^. We hope that future 
studies searching for a divergent non-linear capacitance 
and/or critical spin-spin correlation functions will even- 
tually settle the question about the existence of a ther- 
modynamic glass transition in the 3D lattice model. 

2. Density of states 

The local density of states and the Efros-Shklovskii 
pseudogap has been subject of rather extensive studies in 
the numerical literature^^*^. It has already been shown 
in Ref. 39 that the temperature-dependent distribution 



of instantaneous fields, Pih), reproduces remarkably well 
the numerical results in the high temperature phase of 
the moderately disordered model studied in Ref.[82- Here 
we briefly discuss the low temperature and low energy 
results. 

As we have argued above the fixed point predic- 
tion (|112p for asymptotically low temperatures applies to 
very low energies which are far from reach in numerical 
simulations. However, our results obtained at intermedi- 
ate energies 0.01 < E < 1 compare reasonably well with 
the numerical data of Ref. cf.. Fig. [TT] One should 
however bear in mind that the latter cover only small 
to moderate disorder strengths for which our mean field 
theory is not easy to justify. Nevertheless, as mentioned 
before, the apparent larger pseudogap power a « 2.2—2.4 
in the numerical data is well reproduced by mean field 
theory as a consequence of lattice effects. 

While our theory predicts an essentially disorder in- 
dependent density of states already at moderately low 
energies, the one found in the simulations increases with 
bare disorder. This might be either due to the increas- 
ing difficulty to anneal a sample with stronger disorder 
in a simulation, or it could be a further effect of the lat- 
tice discreteness which affects the non-universal regime 
of energies above T*. 

The theoretical prediction extracted from our lowest 
temperature data for the density of states, P{h; T = 
0.01), lies somewhat below the numerical data. This is 
likely to be ascribed to the clustering of sites with ener- 
gies close to the chemical potential, which was observed 
long ago in simulations of Davies et ali^. This effect is not 
captured correctly in our mean field description (and it is 
also neglected in the ES stability argument). Due to the 
clustering, the numerical data suggests a greater number 
of active sites than there really are, e.g., for the purpose 
of hopping conductivity. Indeed, it is not possible to fill 
two sites in the same small cluster of empty low energy 
sites without paying a lot of energy. We believe that in 
the mean field theory such small clusters effectively occur 
as a single soft site which leads to a lower P{h) than seen 
numerically, while the overall behavior as a function of 
energy is well reproduced. 

However, such a clustering phenomenon should not 
affect the ratio of the density of states at different en- 
ergies and temperatures. In the literature, the ratio 
r = p{E = 0, r = t)/p{E = t,T = 0) was studied numer- 
ically for small t > 0.05 (in energy units of / nt), and 
a rather large value r « 11 was obtained^l by extrapola- 
tion to t = 0. While this ratio could not be reproduced 
by finite T arguments along Efros' self-consistency argu- 
ments^, the present mean field theory comes much closer 
to the numerical finding [with our notation and units 
one has r = P(/i = 0, T = At) /P{h = 2t,T = 0)] . In par- 
ticular our finite temperature results predict the ratios 
r{t = 0.018) ^ 10.8 and r{t = 0.0125) « 8.5, which even- 
tually decreases towards the extreme low energy value 
r{t 0) fa 7.3 which can be extracted from the fixed 
point theory. 
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FIG. 17; (Color online) Symbol lines: 3D data taken from 
Ref. 68, corresponding to uniformly distributed disorder with 
variance W = ((e/2)^) = 0.29,0.58,1.15,2.30 in units of 
e^/4K£ (from bottom to top at low energies). The full lines 
(blue and red) are our mean field predictions for P(/i) at 
T = 0.01, for Gaussian disorder with IF' = 2 and W = 10, 
respectively. Dashed line: Prediction from the ES stability ar- 
gumenlJi, p{E) = (3/7r)£;^. In the plot, we use the standard 
units / K,l and (. for energy and length, respectively. 

3. Dynamics 

It would be interesting to investigate the dynamic as- 
pects of the density of states, or rather, of the distribu- 
tion of time-averaged site magnetizations m = tanh(/3y), 
in order to check for the dynamical self-similarity pre- 
dicted by the replica theory and its interpretation a la 
Sompolinsky. The same applies to the aging dynamics 
in the SK model, where the prediction of dynamical self- 
similarity might serve as a new check for the validity of 
the dynamic interpretation of the mean field solution. 

XI. CONCLUSION 

We have presented a detailed study of the mean field 
theory for electron glasses in 3D. A glass transition is pre- 
dicted once the temperature becomes low enough for fluc- 
tuations in the Thomas-Fermi screening to become criti- 
cally large. The low temperature phase is expected to be 
glassy and remain critical throughout due to a permanent 
marginal stability of typical metastable states. The lat- 
ter also assures the saturation of Efros-Shklovskii's bound 
for the single particle density of states. 

Solving the mean field equations in the low tempera- 
ture phase and in the asymptotic T — > limit, we have 
obtained new quantitative predictions for the evolution 
of the density of states and the compressibility. We have 
given a physical interpretation of the asymptotic fixed 
point occurring in Parisi's flow equations in terms of a 
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dynamical self-similarity of the distribution of time aver- 
aged magnetizations. 

The low T solution exhibits a plateau threshold xp 
in the overlap function Q{x), which we have argued to 
describe the energy landscape in the vicinity of a given 
metastable state. Assuming a dynamics like in the trap 
model, we have used this information to predict the 
asymptotics of the aging function after an excitation. 

This mean field approach is in principle amenable to 
include quantum tunneling between localized states^'^'^. 
The resulting quantum glass will be an interesting candi- 
date to study effects of quantum giassiness on transport 
and slow dynamics. 
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A. L. Efros, L. B. loffe, M. Mezard, Z. Ovadyahu, D. 
Sherrington, and B.I. Shklovskii. This work was in part 
supported by grant PA002-113151 of the Swiss National 
Science Foundation. 



APPENDIX A: MAPPING TO SINGLE SITE 
MODEL 

In this appendix we show how to derive an effective 
local theory using a cumulant expansion. 

Let us consider a Hamiltonian H{(j)) describing a 
generic problem with two body interactions. The fields 
(f) represent the available degrees of freedom (Ising spins 
in our case). We separate the Hamiltonian into a local 
part, Hi{(j)), and a non-local part H2{4>) describing the 
two body interactions 



The variable r is usually a space index, but can also in- 
corporate other coordinates such as a replica label. 

We now perform a cumulant expansion of the partition 
function in terms of the nonlocal term H2{4'), expanding 
around the local limit described by Hi {(f)) . In a graphic 
representation, each diagrams consists of lines represent- 
ing the nonlocal matrix J^{r, r') and vertices. An n-point 
vertex stands for a bare cumulant Af^, defined as 



In (exp {h4)))Hi- 



(A2) 



The diagrammatic expansion is organized in a transpar- 
ent way by introducing the "one-particle irreducible" n- 
point vertex M„(ri, r„) where the denote the space 
(and possibly further internal) indices at which the n ex- 
ternal legs are attached. The two-point correlator can 
then be expressed conveniently as a geometric series. 



G(r,/) = (0(r),^(r')) = + PJ)] 



(A3) 



The effective local theory used in the paper can be 
justified and derived in the limit of large disorder, or al- 
ternatively by invoking a large space dimensionality and 



performing an expansion in 1/D. Analyzing the struc- 
ture of the leading irreducible vertex contributions, one 
finds that their fc-dependence is negligible on the scales 
of interest, i.e., k < Tc- This suggests to treat all irre- 
ducible vertices as local, M„(ri, ...,r„) = unless all r.i 
refer to the same site. In particular, for the local corre- 
lator Eq. (IA3|) simplifies to 



G = G{r,r) 



1 



(A4) 



where we performed a decomposition into the J^q eigen- 
modes q of the matrix J(rr'\ The correlation function 
can simply be written as G = A/2 — M2/S.M2^ defining 
the dressed bond operator 



(A5) 



By restricting to local subdiagrams, the irreducible ver- 
tices can be expanded into a perturbation series in terms 
of the local dressed bond A and the bare cumulants 
{AfO}, Af„=A^„[A,{A^O}]. 

Repeating a similar analysis for a local model 
with Hamiltonian iJ (1°==) (</>) = i7i(0) + H^^°''\<j)), 
we find expressions very similar to the above ones, 
with spatial indices suppressed. In particular, for 



two body interaction H, 



(loc 



find Atf°') 



M. 



(loc) 



[/3A(i°'= 



M 

-1 



[A('-),{A//0}] 



0A(i°'=)0/2, we 
and A(i°"' ^ 



, with the same bare cumulants 

{Af°} and the same functionals A^„ as for the lattice 
model. 

From these observations it becomes clear that the 
lattice correlation functions can be computed with the 
help of the local model on which we impose the self- 
consistency condition A(r, r) = A'^'°'^\ or equivalently, 
M2{r,r) = M^^°''\ From Eq. and its analog in the 

/3A('°'=), we can 



local model, [G^^°'^]-^ = [a4°'Y' 
rewrite the self-consistency equation as 

G = ^(g-1-/3A(i°=)+/3J, 



(A6) 



where G = G(r, r) = G'-^°''\ 

In the main part of the paper we have adopted the 
notations: 



G ^ Q, 


(A7) 


A(i°<=) ^ -(3A, 


(A8) 


a.b 


(A9) 


H2 ^ J2 ''i'^^^'v 


(AlO) 


i>j,a,b 




Hi'""^ ^ -i/3^s»A"''s^ 


(All) 



a.b 
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Further, the irreducible two-point vertex is identified as 
the irreducible polarizability divided by (3, 

M2,ab ^ Tlab/P. (A12) 

In Ref.mthe notation n//3 = (P'^W'^I+'S)-^ was used. 
This was suggestive to indicate two contributions to the 
"self-energy" (3/Tl occurring in the field theory for fields 
conjugate to the spins: a bare quadratic term (3^W^T 
and an additional dressing self-energy S. The latter can 
conveniently be shown to be predominantly local for large 
W. 



APPENDIX B: COMPUTATION OF A ON A 
HYPERCUBIC LATTICE 

In order to compute the Fourier transform of 1 /r^ on 
a hypercubic lattice, we use Ewald type summation 

where 



I{ti,t2;x) 



1 



'dt. 



(B2) 



To obtain well behaved integrals we have to split 
/(O, -l-oo; x) into two terms /(O, e; x) and /(e, -l-oo; x) and 
switch to a momentum sum in the first term. To this end 
we use the representation 

/ df^x 5{x -xi)= [ d'^xY^ [e'^"' - Six)] , (B3) 

where in the right hand side the summation is over the 
vectors of the reciprocal lattice. We then integrate out 
X and change the variable of integration in / to i ^ l/t. 
In this way we find 

— j > (P d-3 



Jk 



G, 



4e 



+ (7) ^ E'^-^w-^V'^" 

where are Misra functions, defined as 

f+00 

0„(z) = / re~''dt. 



(B4) 



(B5) 



Note that for a hypercubic lattice with lattice spacing 
1, one has Gi = 2i:xi. Setting e = tt, we can rewrite 
Eq. jHH as 



Jk — (t)d-3 
2 

+ E 



(B6) 
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In the limit fc — > the first term is dominant and yields 
47r/fc^ in £> = 3 dimensions. 



APPENDIX C: SK MODEL WITH RANDOM 
FIELDS 

The Sherrington Kirkpatrick model is the mean field 
version of spin glasses, describing N spins, all coupled 
mutually with random Gaussian interactions Jij with 
zero mean and variance (J?) — J'^ /N , 



H = 



n SiJijSj 



(CI) 



We also include Gaussian random fields hi with variance 
(/if) = W^, to obtain a model that closely resembles the 
electron glass problem. 

For the SK model the locator approximation"^^ is ex- 
act in the large N limit. The irreducible polarizability 
can easily be shown to be diagonal to leading order in 
N. This allows for the same local approximation as we 
have used in the electron glass problem, the only differ- 
ence being that the interaction matrix has a different 
eigenvalue spectrum. 



p{X)dX = —y/J^ 
ttJ 



A2e(J- |A|)dA, 



(C2) 



which is normalized such that J p{X) dX ~ I. 

This replaces the " mode density" of the electron glass 



p(A) dX = 



d^k 



6{X - Jk) dX. 



By mapping to a single site model of the form 

^^o = -yE^"(^«'' + ^')^'' 

ab 



(C3) 



(C4) 



one obtains the self-consistency condition for the SK 
model with random fields. 



Qab 



{SaSb/Hn 



dXp{X) 



f3X' 
(3^ A 



ab 



(C5) 



In this case, the integral over A can be done explicitly. 
After a little algebra one recovers the saddle point equa- 
tion from the standard replica treatment which imposes 
that the overlap matrix and the effective coupling matrix 
coincide, Qab = ^ab- 

The further steps to solve (|C5[) are analogous to the 
case of the electron glass. In particular, the glass transi- 
tion is given by the solution of 



f3' 



cosh (Py) 



(C6) 



with 



PRsiy; Qrs) 



Qns = J dyPRs(y;QRs)tanh2(/3y), (C7) 
exp[-yV2(iy2_^QRs)] 



(08) 
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FIG. 18: (Color online) Universality of the SK model in ran- 
dom fields. Top: Scaling of the local field distribution P{y) for 
various temperatures 10 < /3 < 60 and disorders < W < 1.5. 
Bottom: Unsealed distribution P{y) at /3 = 20 for various dis- 
order strengths, W = 0, 0.3, 0.6, 0.9, 1.2. 



From these equations it is easy to establish the scaling of 
the Almeida Thouless-line Tc ^ 1/W a.t large disorder. 

In Fig. [18] we illustrate the universality (disorder- 
independence) of the linear pseudogap, and its scaling 
behavior P{y) = T^sv.[y/T) at low temperature. 



APPENDIX D: DERIVATION OF THE ONSAGER 
TERM 

The Onsager term ho represents the response of the 
environment to a spin at a given site 0. The latter should 
thus be excluded from all ring diagrams in ((23)) . In order 
to establish a relation between ho and the self-consistent 
coupling matrix Kab we analyze the diagrammatic expan- 
sion for the local overlap 

PQ^:^ = [n-n^n + njn^n- + ...]^, (di) 

where the matrix products in the brackets involve sums 
over replicas and all sites, including i. Hah = f3M2.ab 
is the irreducible local polarizability or two point vertex. 
The result is independent of the site i due to the disorder 
average. 

Let us now define the "Onsager matrix" 



where all internal polarizabilities are restricted to sites 
j ^ i- With this matrix we can rearrange the sum 
(|D1[) according to the number of intermediate visits of 
the site i. 



PQ'^^ = n + Wion + nhoUhoU 
1 



n-1 - hr 



(D3) 



ab 



Comparing to Eq. ([5]) we see that in the mean field ap- 
proximation, the Onsager matrix is essentially just the 
couphng matrix, ho,ab = P^ab- 

In order to obtain the (average) Onsager back reaction 
in a given metastable state, we observe that the local, 
connected spin correlator calculated within a single state 
is given by Qaa - Qa^b = Q - Q{x = 1) = [Q]{x = 1). 
The associated Onsager term is therefore 



ho = ho,aa-ho,a^b = H^)] 



= Xi ^{i^c) 



2y/TrKc, 



(D4) 



where we have used Eqs. (|57l60p in the second line and 
the asymptotic expression ([9]) for the last step. 

Note that for the SK model where A = Q one correctly 
recovers ho = f3[Q - Q(l)] - /?(! - Qea). 



APPENDIX E: REPLICA FOURIER 
TRANSFORM 

The replica Fourier transform***'^^ performs the diago- 
nalization of a generic Parisi matrix Q to a representation 
in which it is straightforward to compute any analytic 
function of the matrix. 

An ultrametric matrix with K steps of replica sym- 
metry breaking is fully described by its diagonal en- 
tries Qaa = Q and off-diagonal blocks with entries 
Qk, Qk-1, ■ ■ ■ , Qi, Qq in order of their distance from the 
diagonal. Such a matrix can be written as 

K 

Q - QoUo + Y,{Qk- Qk-i)n,n, + {Q~ Qm^ ) id, (El) 
fc=i 

where < mi < • • • < rriK < mn+i = 1 are the sizes 
of the replica clusters in the hierarchy, and TZm denotes 
a n X 71 matrix with n/m blocks on the diagonal each of 
which has all entries equal to one. 

To compute functions of such a matrix, it is convenient 
to introduce the replica Fourier transform 



K 

El 

fc=l 



. Ak + Qc id, 



(E2) 



with the matrices Ak = TZmt+i/nT-k+i —Ti-mk/'n^k for k 
1, . . . , K, which satisfy orthonormality relations 



[juj - jnjnj 



lb ' 



(D2) 



(E3) 



27 



This makes it is easy to compute any analytic function 
G = f{Q) by Taylor expansion. 

In the representation (|E2|1 one easily finds 



S = f'{Qc)Qono 



(E4) 



K 



The relation between [Q]f. and Qk is 



= Qo + E[Q],f^ —) + ^^, (E5) 



K 



implying 



Wlk TOfe+i 



Qk ~ Qc — [Q]k ■ 



(E7) 
(E8) 



In the continuum limit [K — + cx3, Qfe — ^ Q[x — rn/c)) 
Eqs. (113) turn into 



d[Q] _ dQ 
dx dx 



(E9) 



Integrating, we find the continuous replica Fourier trans- 
form 

[Q]{x) = xQ{x)- f Q(y)dy, (ElO) 
Jo 

Qc ^ Q- f Q{x)dx. (Ell) 
Jo 

The continuous Fourier components oi Q = f{Q) can 
be read off from (IE4D as 



Gc = f{Qc), (E12) 
[G] (x) = fiQc) - f{Qc - [Q] (x)), (E13) 
Go = f'{Qc)Qo, (E14) 

which implies the useful identity 

- [G] (a:) = /(Q, - [Q] (a:)), (E15) 

which is used repeatedly in the main part of the paper. 



APPENDIX F: NUMERICAL 
IMPLEMENTATION OF PARISI'S EQUATIONS 

In this appendix we describe several details of the nu- 
merical integration of the flow equations (|49l5ip . which 



are subtle in particular in the low temperature limit. The 
integration of the flow equations forms the crucial part 
in the iterative scheme outlined in the main text. 

Given A(x) we first perform the integration of the flow 
for m{x,y). Since the significant features of m live on 
the scale y l/{[3x), it is convenient to transform to the 
new variable z = f{l3x)y with f{b) ^ b for large b. We 
have used f(b) = b + c with c a constant of order 10, 
which yields a smooth behavior of m(x, z), also at low x. 

In a second step we integrate the flow of the distribu- 
tion function P{x,y). This is tricky at low temperature 
because there are at least two disparate energy scales, 

I kI and T, which are both characteristic scales for 
P(j/). If the disorder W is strong, it introduces an even 
larger scale, on which varies significantly. 

This problem can be solved by noting that to high 
accuracy we can set mix^ y) — sgn(?/) when f{l3x)\y\ ^ 1, 
where m approaches ±1 exponentially. In this region of 
fields the flow equation (|5ip simplifies to the linear partial 
differential equation 



P"{y;x) 



P{y;x) = A{x) 



where s = sgTi{y). Its kernel is 



3 l3xP\y; x) 



(Fl) 



G{x,y;x',y') = 



exp 



iy-y'-sf3 f^, uA{u)du)^ 



2[A{x}-A{x')] 



^27T[Aix) ~ A(x')] 



(F2) 



with the help of which we can propagate initial conditions 
from P{x = 0, y) to finite x in order to obtain boundary 
conditions for large y ^ 1/ f{(3x). A closer inspection 
shows, that in order to provide boundary conditions for 
all x, one has to be a bit more careful so as not to leave 
the range of validity of the linear Eq. (jFip . i.e., the range 
where |m| > 1 — e for some fixed e (we usually imposed at 
least e < 10""). The latter requires a stepwise forward- 
integration from = to xi, from xi to X2 etc., where 
in each interval [xn,Xn+i\ the boundary conditions are 
imposed at a different (decreasing) ?/„. 

It is usually sufficient to store P(x„, \y\ < 2?/„) both to 
calculate boundary conditions for P(x„ < x < Xn^i,yn) 
and new initial conditions P{xn+i,\y\ < 2yn+i). The 
partitioning {Xn} is determined dynamically in such a 
way as to ensure that \m{x,yn — A„)| > 1 — e for all 
Xn < X < Xn+i which is required to justify the use of 
Eq. jFT])- Here, A„ = /3 duuk(u) 

With this technique, the whole interval < a; < 1 can 
be covered maintaining a high precision in the regime of 
interest f{(3x)y = 0(1) for the purpose of calculating the 
overlap Q{x) and to impose self-consistency. We could 
reach as low temperatures as /? = 100 without encounter- 
ing numerical problems. In this way we have achieved a 
high precision for the universal tails in A(5), Q(6), rsc(6) 
for 6 > 1. 
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